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Oh' Abstract High-energy X-rays and y-rays from solar flares were discovered just over fifty 



years ago. Since that time, the standard for the interpretation of spatially integrated flare 
X-ray spectra at energies above several tens of keV has been the collisional thick-target 
model. After the launch of the Reuven Ramaty High Energy Solar Spectroscopic Imager 
(RHESSI) in early 2002, X-ray spectra and images have been of sufficient quality to allow a 
greater focus on the energetic electrons responsible for the X-ray emission, including their 
origin and their interactions with the flare plasma and magnetic field. The result has been 
new insights into the flaring process, as well as more quantitative models for both electron 
acceleration and propagation, and for the flare environment with which the electrons interact. 
In this article we review our current understanding of electron acceleration, energy loss, and 
propagation in flares. Implications of these new results for the collisional thick-target model, 
Q ■ for general flare models, and for future flare studies are discussed. 

H ' 

Keywords Sun: flares, Sun: X-rays, gamma rays, Sun: radio radiation 



Oh- 



> 



Contents 



1 Introduction 

2 Thin- and thick-target X-ray emission 

3 Low-energy cutoffs and the energy in non-thermal electrons 

3.1 Why do we need to determine the low-energy cutoff of non-thermal electron distributions? . 

3.2 Why is the low-energy cutoff difficult to determine? 

• 3.3 What is the shape of the low-energy cutoff, and how does it impact the photon spectrum and 

PmhP- [2 



1 NASA Goddard Space Flight Center, Code 671, Greenbelt, MD 20771, U.S.A. E-mail: Gor- 
don. D.Holman@nasa.gov 

2 Lockheed Martin Advanced Technology Center, Solar and Astrophysics Laboratory, Organization ADBS, 
Building 252, 3251 Hanover Street, Palo Alto, CA 94304, U.S.A. 

■ 3 Astrophysikalisches Institut Potsdam 

J.^ ' 4 Department of Physics & Astronomy, University of Glasgow, Glasgow, G12 8QQ, Scotland, UK 

5 Harvard-Smithsonian Center for Astrophysics P-148, 60 Garden St., Cambridge, MA 02138, U.S.A. 

6 Department of Physics and Astronomy, University of Glasgow, Kelvin Building, Glasgow, G12 8QQ, U.K. 

7 Space Sciences Lab, UC Berkeley, CA, U.S.A. 

8 Department of Computing and Mathematics, University of Bradford, Bradford, BD7 1DP, UK 



2 



Holman et al. 



3.4 Important caveats 

3.5 Determinations of E c and electron energy content from flare data 

4 Nonuniform ionization in the thick-target region 

4.1 Electron energy losses and X-ray emission in a nonuniformly ionized plasma 

4.2 Application to flare X-ray spectra 

5 Return current losses 

5.1 The return current electric field 

5.2 Impact on hard X-ray spectra 

5.3 Observational evidence for the presence of the return current 

6 Beam-plasma and current instabilities 

7 Height dependence and size of X-ray sources with energy and time 

7.1 Footpoint Sources 

7.2 Loop Sources and their Evolution 

8 Hard X-ray timing 

8.1 Time-of-Flight Delays 

8.2 Trapping Delays 

8.3 Thermal Delays 

8.4 Multi-Thermal Delay Modeling with RHESSI 

9 Hard X-ray spectral evolution in flares 

9.1 Observations of spectral evolution 

9.2 Interpretation of spectral evolution 4(] 

10 The connection between footpoint and coronal hard X-ray sources 

10.1 RHESSI imaging spectroscopy 

10.2 Relation between coronal and footpoint sources 

10.2.1 Observed difference between coronal and footpoint spectral indices 

10.2.2 Differences between footpoints 

10.3 Spectral evolution in coronal sources 

10.4 Interpretation of the connection between footpoints and the coronal source 

1 1 Identification of electron acceleration sites from radio observations 

12 Discussion and Conclusions 

12.1 Implications of X-ray observations for the collisional thick-target model 

12.2 Implications of X-ray observations for electron acceleration mechanisms and flare models 

12.3 Implications of current results for future flare studies in hard X-rays 



1 Introduction 

A primary characteristic of solar flares is the acceleration of electrons to high, suprathermal 
energies. These electrons are observed directly in interplanetary space, and indirectly a t 
the Sun through the X-ray, y-ray, and radio emissions they emit ( iHudson & Rvaniri995h . 
Understanding how these electrons are produced and how they evolve is fundamental to 
obtaining an understanding of energy release in flares. Therefore, one of the principal goals 
of solar flare research is to determine when, where, and how these electrons are accelerated 
to suprathermal energies, and what happens to them after they are accelerated to these high 
energies. 

A major challenge to obtaining an understanding of electron acceleration in flares is 
that the location where they are accelerated is not necessarily where they are most easily 
observed. The flare-accelerated electrons that escape the Sun are not directly observed until 
they reach the instruments in space capable of detecting them, usually located at the distance 
of the Earth from the Sun. The properties of these electrons are easily modified during thei r 
long journey from the flaring region to the detecting instruments (e.g., lAgueda et al.l E009). 
Distinguishing flare-ac celerated elec trons from electrons accelerated in interplanetary shock 
waves is also difficult (Kahler 2007). 

The electrons that are observed at the Sun through their X-ray or y-ray emissions radiate 
most intensely where the density of the ambient plasma is highest (see Section^. Therefore, 
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the radiation from electrons in and near the acceleration region may not be intense enough to 
be observable. Although these radiating electrons are much closer to the acceleration region 
than those detected in interplanetary space, their properties can still be significantly modified 
as they propagate to the denser regions where they are observed. The radio emission from 
the accelerated electrons also depends on the plasma environment, es pecially the magneti c 
field strength for the gyrosynchrotron radiation observed from flares dBastian et al.| [l998). 
Therefore, determining when, where, and how the electrons were accelerated requires a 
substantial amount of deductive reasoning. 

Here we focus primarily on the X-ray emission from the accelerated electrons . Interplan- 
etary electrons and low-energ y emissions are add ressed bv lFletcher et alJ J201 lb. whi le the 
y-ray emission is addressed by Vilm er et al.l d201ll) . and the radio by W hite et al. 1 201 lh - The 
X-rays are predominantly electron-ion bremsstrahlung (free-free radiation), emitted when 
the accelerated electrons scatter off ions in the ambient thermal plasma. Issues related to the 
emission mechanism and deducing the properties of t he em itting electrons from the X-ray 
observations are primarily addressed bv lKontar et alJ ( 1201 ll) . Here we address the interpre- 
tation of the X-ray observations in terms of flare models, and consider the implications of 
the observations for the acceleration process, energy release in flares, and electron propaga- 
tion. Specific models for particle acceleration and energy release in flares are addressed by 
IZharkovaetaflfcOllh . 

The accelerated electrons interact with both ambient electrons and ions, but lose most of 
their energy through electron-electron Coulomb collisions. Consequently, the brightest X- 
ray sources are associated with high collisional energy losses. These losses in turn change 
the energy distribution of the radiating electrons. When the accelerated electrons lose their 
suprathermal energy to the ambient plasma as they radiate, the source region is called a 
thick target. Electrons streaming downward into the higher densities in lower regions of the 
solar atmosphere, or trapped long enough in lower density regions, will emit thick-target 
X-rays. Hence, thick- target models are important to understanding the origin and evolution 
of accelerated electrons in flares. Thick-target X-ray emission is addressed in Section|2] 

The total energy carried by accelerated electrons is important to assessing acceleration 
models, especially considering that these electrons carry a significant fraction of the total 
energy released in flares. Also, the energy carried by electrons that escape the acceleration 
region is deposited elsewhere, primarily to heat the plasma in the thick-target source regions. 
The X-ray flux from flares falls off rapidly with increasing photon energy, indicating that the 
number of radiating electrons decreases rapidly with increasing electron energy. Therefore, 
the energy carried by the accelerated electrons is sensitive to the value of the low-energy 
cutoff to the electron distribution. The determination of this low-energy cutoff and the total 
energy in the accelerated electrons is addressed in Section[3] 

In the standard thick-target model, the target plasma is assumed to be fully ionized. 
If the target ionization is not uniform, so that the accelerated electrons stream down to 
cooler plasma that is partially ionized or un-ionized, the X-ray spectrum is modified. This is 
addressed in Section[4] 

Observations of the radiation from hot flare plasma have show n this plasma to prim arily 
be confined to magnetic loops or arcades of magnetic loops (cf. lAschwandenl 2004). The 
observations also indicate that the heating of this plasma and particle acceleration initially 
occur in the corona above these hot loops (see Section [T2.2l and lFletcher et al.ll201 lh . When 
the density structure in these loops is typical of active region loops, or at least not highly en- 
hanced above those densities, the highest intensity, thick-target X-ray emission will be from 
the footpoints of the loops, as is most often observed to be the case. If accelerated elec- 
trons alone, unaccompanied by neutralizing ions, stream down the legs of the loop from the 
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acceleration region to the footpoints, they will drive a co-spatial return current in the ambi- 
ent plasma to neutralize the high current associated with the downward-streaming electrons. 
We refer to this primarily downward-streaming distribution of energetic (suprathermal) elec- 
trons as an electron beam. The electric field associated with the return current will decelerate 
electrons in the beam, which can in turn modify the X-ray spectrum from the accelerated 
electrons. This is addressed in Section[5] 

Both the primary beam of accelerated electrons and the return current can become unsta- 
ble and drive the growth of waves in the ambient plasma. These waves can, in turn, interact 
with the electron beam and return current, altering the energy and angular distributions of 
the energetic electrons. These plasma instabilities are discussed in Section|6] 

The collisional energy loss rate is greater for lower energy electrons. Therefore, for 
suprathermal electrons streaming downward to the footpoints of a loop, the footpoint X-ray 
sources observed at lower energies should be at a higher altitude than footpoint sources ob- 
served at higher X-ray energies. The height dispersion of these sources provides information 
about the height distribution of the plasma density in the footpoints. The spatial resolution 
of the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI — see Lin et al. 
2002) has made such a study possible. This is described in Section|7] RHESSI has observed 
X-ray sources that move downward from the loop top and then upward from the footpoints 
during some flares. This source evolution in time is also discussed in Section|7] 

If electrons of all energies are simultaneously injected, the footpoint X-ray emission 
from the slower, lower energy electrons should appear after that from faster, higher energy 
electrons. The length of this time delay provides an important test for the height of the accel- 
eration region. Longer time delays can result from magnetic trapping of the electrons. The 
evolution of the thermal plasma in flares can also exhibit time delays associated with the bal- 
ance between heating and cooling processes. These various time delays and the information 
they provide are addressed in Section[8] 

An important diagnostic of electron acceleration and propagation in flares is the time 
evolution of the X-ray spectrum during flares. In most flares, the X-ray spectrum becomes 
harder (flatter, smaller spectral index) and then softer (steeper, larger spectral index) as the 
X-ray flux evolves from low to high intensity and then back to low intensity. There are 
notable exceptions to this pattern, however. Spectral evolution is addressed in Section[9] 

One of the most important results from the Yohkoh mission is the discovery in some 
flares of a hard (high energy) X-ray source above the top of of the thermal (low energy) 
X-ray loops. This, together with the Yohkoh observations of cusps at the top of flare X- 
ray loops, provided strong evidence that energy release occurs in the corona above the hot 
X-ray loops (for some flares, at least). Although several models have been proposed, the 
origin of these "above-the-looptop" hard X-ray sources is not well understood. We need to 
determine how their properties and evolution compare to the more intense footpoint hard 
X-ray sources. These issues are addressed in SectionfTOl 

Radio observations provide another view of accelerated electrons and related flare phe- 
n omena. Although ra dio emission and its relationship to flare X-rays are primarily addressed 
in White et al.l feOllh . some intriguing radio observations that bear upon electron accelera- 
tion in flares are presented in Section [TT] 

RHESSI observations of flare X-ray emission have led to substantial progress, but many 
questions remain unanswered. Part of the progress is that the questions are different from 
those that were asked less than a decade ago. The primary context for interpreting the X-ray 
emission from suprathermal electrons is still the thick-target model, but the ultimate goal is 
to understand how the electrons are accelerated. In Section[l2]we summarize and discuss the 
implications of the X-ray observations for the thick-target model and electron acceleration 
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mechanisms, and highlight some of the questions that remain to be answered. Implications 
of these questions for future flare studies are discussed. 



2 Thin- and thick-target X-ray emission 

As was summarized in Section [TJ the electron-ion bremsstrahlung X-rays from a beam of 
accelerated electrons will be most intense where the density of target ions is highest, as 
well as where the flux of accelerated electrons is high. The local emission (emissivity) at 
position r of photons of energy e by electrons of energy E is given by the plasma ion 
density, n(r) ions cm -3 , times the electron-beam flux density distribution, F(E,r) electrons 
cm -2 s _1 keV -1 , times the differential electron-ion bremsstrahlung cross-section, Q(e,E) 
cm 2 keV _I . For simplicity, we do not consider here t he angular distribut ion of the beam 
electrons or of the emitted photons, topics addressed in Kontar et al.l d201 lh . 

The emissivity of the radiation at energy e from all the electrons in the beam is obtained 
by integrating over all contributing electron energies, which is all electron energies above 
the photon energy. The photon flux emitted per unit energy is obtained by integrating over 
the emitting source volume (V) or, for an imaged source, along the line of sight through 
the source region. Finally, assuming isotropic emission, the observed spatially integrated 
flux density of photons of energy e at the X-ray detector, 7(e) photons cm~ 2 s _1 keV -1 , is 
simply the flux divided by the geometrical dilution factor AnR 2 , where R is the distance to 
the X-ray detector. Thus, 

/ ( £ ) = 4^2 / f ' n{r)F{E,r)Q{e,E)dEdV. (2.1) 

We refer to 7(e) as the X-ray flux spectrum, or simply the X-ray spectrum. The spectrum 
obtained directly from an X-ray detector is generally a spectrum of counts versus energy 
loss in the detector, which must be converted to an X-ray flux spect rum by correcting for the 
detector response as a function of photon energy (see, for example. ISmith et aDl2002h . 

Besides increasing the X-ray emission, a high plasma density also means increased 
Coulomb energy losses for the beam electrons. In a plasma, the bremsstrahlung losses are 
small compared to the collisional losses to the plasma electrons. For a fully ionized plasma 
and beam electron speeds much greater than the mean speed of the thermal electrons, the 
(nonrelativistic) energy loss rate is 

dE/dt = -{K/E)n e (r)v{E), (2.2) 

where K = 2ne 4 A ee , A ee is the Coulomb logarithm for electron-electron collisions, e is the 
electron charg e, n e (r) is the plasma elec tron number density, and v[E) is the speed of the 



electron (see Brown|[l97 lb lEmsliel lT978"). The coefficient K is usually taken to be constant, 



although A ee depends weakly on the electron energy and plasma density or magnetic field 
strength, typically falling in the range 20 - 30 for X-ray-emitting electrons. Taking a value 
of 23 for A ee , the energy loss rate in keV s _1 or erg s _1 with E in keV is numerically 
determined by 

K = 3.00 x 10~ 18 \ keVcm 2 = 4.80 x 10~ 27 J ergcm 2 . (2.3) 

Here and in equations to follow, notation such as (A cc /23) is used to show the scaling of 
a computed constant (here, K) with an independent variable and the numerical value taken 
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for the independent variable (A ee = 23 in Equation 12.3b . If the plasma is not fully ionized, 
K also depends on the ionization state (see Section l4~Tt . 

Noting that vdt = dz, Equation l2.2l can be simplified to dE/dN e = —K/E, where dN e (z) = 
n e (z)dz and N e (z) (cm -2 ) is the plasma electron column density. (Here we treat this as a 
one-dimensional system and do not distinguish between the total electron velocity and the 
velocity component parallel to the magnetic field.) Hence, the evolution of an electron's 
energy with column density is simply 

E 2 =E^-2KN e , (2.4) 

where Eq is the initial energy of the electron where it is injected into the target region. 
For example, a 1 keV electron loses all of its energy over a column density of 1/ (2K) = 
1.7 x 10 17 cm~ 2 (for A ee = 23). A 25 keV electron loses 20% of its energy over a column 
density of 3.8 x 10 19 cm~ 2 . 

If energy losses are not significant within a spatially unresolved X-ray source region, 
the emission is called thin-target. If, on the other hand, the non-thermal electrons lose all 
their suprathermal energy within the spatially unresolved source during the observational 
integration time, the emission is called thick-target. We call a model that assumes these 
energy losses are from Coulomb collisions (equations 12.21 & [2731 ) a collisional thick-target 
model. Collisional thick-target models h ave been applied to flare x - ray/y -ray emission since 
the discovery of this emission in 1958 dPeterson & WinckleJ [l958lll959h . 

The maximum information that can be obtained about the accelerated electrons from an 
X-ray spectrum alone is contained in the mean electron flux d istribution, the plasma-density- 
weigh ted, target-averaged electron flux density distribution terown et al.l2003l : lKontar et alJ 
l201ll) . The mean electron flux distribution is defined as 

F(E) = — [ n(r)F(E,r)dV electrons cm~ 2 s _1 keV -1 , (2.5) 
nV Jy 

where n and V are the mean plasma density and volume of the emitting region. As can be 
seen from Equation 12. II the product hVF can, in principle, be deduced with only a knowl- 
edge of the bremsstrahlung cross-section, Q(e,E). Additional information is required to de- 
termine if the X-ray emission is thin-target, thick-target, or something in between. The flux 
distribution of the emitting electrons and the mean electron flux distribution are equivalent 
for a homogeneous, thin-target source region. 

Equation 12.11 gives the observed X-ray flux in terms of the accelerated electron flux 
density distribution throughout the source. However, we are interested in the electron dis- 
tribution injected into the source region, Fq(Eq,To), since that is the distribution produced 
by the unknown acceleration mechanism, including any modifications during propagation 
to the source region. To obtain this, we need to know how to relate F(E,r) at all locations 
within the source region to Fo(iio,ro). Since we are interested in the X-rays from a spa- 
tially integrated, thick-target source region, the most direct approach is to first compute the 
bremsstrahlu ng photon yie ld from a single electron of energy Eq, v(e,£b) photons keV -1 
per electron ferown|[l97 ll) . As long as the observational integration time is longer than the 
time required for the electrons to radiate all photons of energy e (i.e., longer than the time 
required for energy losses to reduce all electron energies to less than e), the thick-target 
X-ray spectrum is then given by 

I,hick ^ = 4^l ^o(E )v(e,E )dEo, (2.6) 
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where j2"o(£b) is the electron beam flux distribution (electrons s keV ). ^q(Eq) is the 
integral of Fo (•*<)> -So) over the area at the injection site through which the electrons stream 
into the thick-target region. 

The rate at which an electron of energy E radiates bremsstrahlung photons of energy 
e is n(r)v(E)Q(e,E). The photon yield is obtained by integrating this over time. Since the 
electrons are losing energy at the rate dE/dt, the time integration can be replaced by an 
integration over energy from the initial electron energy Eq to the lowest energy capable of 
radiating a photon of energy e: 

v(e,E )=r « r ^%f t ^ dE . (2.7) 
Je dE/dt 



Using Equation I2.2l for dE/dt, Equation l2.6l becomes 

i i r r E ° 



'thick 



(e) = ^ s£. Wo) L EQ ^ dEdE - (2 - 8) 



We have used the relationship n e = Zn, where Z ~ 1.1 is the ion-species-number-density- 
weighted (or, equivalently, relative-ion-abundance-weighted) mean atomic number of the 
target plasma. Thus, the thick- target X-ray flux spectrum does not depend on the plasma 
density. However, the plasma must be dense enough for the emission to be thick-target, 
i.e., dense enough for all the electrons to be thermalized in the observation time interval. 
Integration of Equation l2.2l shows that this typically implies a plasma density > 10 11 — 10 12 
cm -3 for an observational integration time of 1 s (see Sections flO.4l and fl2.H for more about 
this). This condition is well satisfied at loop footpoints. 

Observed non-thermal X-ray spectra from solar flares can usually be well fitted with a 
model photon spectrum that is either a single or a double power-law. For a single power-law 
electron flux distribution of the form ^(E) = AE~ S , the photon spectrum is also well ap- 
proximated by the power-law form 7(e) = Iq£^ 7 . The relationship between the electron and 
photon spectral indices 8 and yean most easily be obtained from equations 12. H and l2.8l us- 
ing the Kramers approximation to the nonrelativistic Bethe-Heitler (NRBH) bremsstrahlung 
cross section (see Koch & Motz 1959 for bremsstrahlung cross-sections). The NRBH 
cross-section is given by: 

Q NmH (e,E) = ^J ] + ^ZlIl ) cm 2 keV-\ (2.9) 
£E \ 1 - a/ 1 - e /E J 



where go = 7.90 x 10 -25 cm 2 keV andZ 2 ~ 1.4 is the ion-species-number-density-weighted 
mean square atomic number of the target plasma. The Kramers approximation to this cross- 
section is Equation 12.91 without the logarithmic term. The bremsstrahlung cross-section is 
zero for e > E, since an electron cannot radiate a photon that is more energetic than the elec- 
tron. Analytic expressions for the photon flux from both a uniform thin-target source and a 
thick-target source can be obtained with the Kramers and the NRBH c ross-sections when the 
electr on flux distribution has the single-power-law form l lBrownl 1 97 lllTandberg -Hanssen & Emsliel 
Il988f) . The thin-target result also generalizes to the photon flux from a single-power-law 
mean electron flux distribution. 

For a uniform thin-target source and ,^(E~) = AE~ S , 

We) = 3 - 93 x 10_52 (-^r) ( Y4 ) NA Pn,(S)£- {5+1 \ (2.10) 
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giving y t ta„ = 5 + 1. The photon energy e is in keV, the distance from the source to the X-ray 
detector is taken to be one Astronomical Unit (1 AU), a typical value of 1.4 is taken for Z 2 , 
and TV is the ion column density. The power-law-index-dependent coefficient for the NRBH 
cross-section is 

5(5,1/2) 

Ptn{0) = ^ , (2.11) 

where B(x,y) is the standard Beta function. In the Kramers approximation, j6 f „(5) = 1/5. 
Typical values for the ion column density and A, the differential electron flux at 1 keV, are 
N = 10 18 -10 20 cm" 2 and A = 10 34 -10 38 electrons s" 1 keV" 1 . 
For a thick-target source region, 

-»i±™\*(*j2\ fa,,.,™, 



We) = 1.17x10-^— ) [j^j^jA^e-^, (2.12) 

giving Yihick = 8 — 1. The power-law-index-dependent coefficient for the NRBH cross-section 
is 

5(5-2,1/2) 

fe(gH (5-l)(5-2) " (2J3) 
In the Kramers approximation, ftt ( 5 ) = l/[(<5 — 1 ) ( ^ — 2) ] . 



CO 
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Fig. 2.1 The power-law-index-dependent terms j8/„(5) (equation 12. Ill and Ptk(S) (equation 12.131 in the 
analytic expressions for the thin-target and the thick-target photon flux from a power-law electron flux dis- 
tribution (equations 12.1 0l and l2. 12) . The solid curves are for the nonrelativistic Bethe-Heitler bremsstrahlung 
cross-section and the dotted curves are for the Kramers cross-section. Note that the j3 axis is linear for the 
thin-target coefficient, and logarithmic for the thick-target coefficient. 



The coefficients ft,, (5) and j6^(5) for both the NRBH and the Kramers cross-sections 
are plotted as a function of 5 in Figure 12.11 The Kramers and NRBH results are equal for 
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thin-target emission when 8 ~ 3.4, and for thick-target emission when 8 ~ 5.4. For the 
plotted range of 8, the Kramers approximation can differ from the NRBH result by over 
90%. For 8 in the range 3-10, the Kramers result can differ from the NRBH result by as 
much as 76% and 57% for thin- and thick-target emission, respectively. 

It is important to recognize that the above power-law relationships are only valid if the 
electron density distribution, F(r,E) electrons citT 2 s _1 keV -1 , or the electron flux dis- 
tribution, JF(E) electrons s _1 keV -1 , is assumed to have a power-law energy dependence. 
It is sometimes convenient to work with the electron density distribution, f(r,E) (electrons 
cm~ 3 keV~'), rather than the flux density distribution, especially when considering thin- 
target emission alone or comparing X-ray spectra with radio spectra. The flux density and 
density distributions are related through F(r,E) = f(r,E)v(E). If the electron density dis- 
tribution rather than the flux or flux density distribution is assumed to have a power-law 
index 8', so that f(r,E) °c E~ s , the relationships between this power-law index and the 
photon spectral index become y,/„>, = 8' + 0.5 and Ytkick = 8' — 1.5. 

The simple power-law relationships are not valid if there is a break or a cutoff in the elec- 
tron distribution at an energy less than ~2 orders of magnitude above the photon energies 
of interest. Since all electrons with energies above a given photon energy e contribute to the 
bremsstrahlung at that photon energy, for the power-law relationships to be valid the break 
energy must be high enough that the deficit (or excess) of electrons above the break energy 
does not significantly affect the photon flux at energy e. The power-law relationship is typ- 
ically not accurate until photon energies one to two orders of magnitude below the break 
energ y, depending on the steepness of the power-law electron distribution (see Figures 9 & 
10 of lHolmarj r2003). Thus, for example, these relationships are not correct for the lower 
power-law index of a double power-law fit to a photon spectrum at photon energies within 
about an order of magnitude below the break energy in the double power-law electron dis- 
tribution. Equation l2.1l or l2~8l can be used to numerically compute the X-ray spectrum from 
an arbitrary flux distribution in electron energy. 

When electrons with kinetic energies approaching or exceeding 511 keV significantly 
contribute to the radiation, the rela tivistic Bethe-Heitler bremsst rahlung cros s-section (Equa- 
tion 3BN of iKoch & Motzlll959l) or a close approximation dHaud 119971) must be used. 
lHaud dl997h has shown that the maximum error in the NRBH cross-section relative to 
the relativistic Bethe-Heitler cross section becomes greater than 10% at electron energies 
of 30 keV and above. Numerical computations using the relativistic Bethe-Heitler cross- 
section have been incorporated into the RHESSI spectral analysis software (OSPEX) for 
both thin- and thick-target emission from, in the most general case, a broken-power-law elec- 
tron flux distribution with both low- and high-energy cutoffs (the functions labeled "thin" 
and "thick" u sing the IDL programs brm_bremspec . pro and brm_bremthick.pro - see 
lHolmanll2003h . Faster versions of these programs are now available in OSPEX and are cur- 
rently labeled "thin2" and "thick2" and use the IDL programs brm2_thintarget .pro and 
brm2_thicktarget . pro. 

The analytic results based on the NRBH cross-se ction have been gen eralized to a broken- 
power-law electron flux distribution with cutoffs by iBrowne t al. (2008). They find a max- 
imum error of 35% relative to results obtained with the relativistic Bethe-Heitler cross- 
section for the range of parameters they consider. These results provide the fastest method 
for obtaining thin- and thick-target fits to X-ray spectra in the RHESSI spectral analysis soft- 
ware, where they are labeled "photon_thin" and "photon_thick" and use the IDL programs 
f _photon_thin.pro and f _photon_thick.pro. 
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3 Low-energy cutoffs and the energy in non-thermal electrons 

One of the most important aspects of the distribution of accelerated electrons is the low- 
energy cutoff. The acceleration of charged particles out of the thermal plasma typically 
involves a competition between the collisions that keeps the particles thermalized and the 
acceleration mechanism. The particles are accelerated out of the tail of the thermal distribu- 
tion, down to the lowest particle energy for which the acceleration mechanism can overcome 
the collisional force. Thus, the value of the low-energy cutoff can provide information about 
the force of the acceleration mechanism. More generally, as discussed below, the electron 
distribution must have a low-energy cutoff (1) so that the number and energy flux of elec- 
trons is finite and reasonable, and (2) because electrons with energies that are not well above 
the thermal energy of the plasma through which they propagate will be rapidly thermalized. 
Knowledge of the low-energy cutoff and its evolution during a flare is critical to determin- 
ing the energy flux and energy in non-thermal electrons and, ultimately, the efficiency of the 
acceleration process. 

3.1 Why do we need to determine the low-energy cutoff of non-thermal electron 
distributions? 

An important feature of the basic thick-target model is that the photon spectrum /(e) is di- 
rectly determined by the injected electron flux distribution ^q(Eq), As can be seen from 
Equation 12.81 no additional parameters such as source density or volume need to be deter- 
mined. Consequently, by integrating over all electron energies, we can also determine the 
total flux of non-thermal electrons, N,,,/, electrons s _1 , the power in non-thermal electrons, 
Pnth er g s -1 , and, integrating over time, the total number of, and energy in, non-thermal 
electrons. 

The total non- thermal electron number flux and power are computed as follows: 

N nlh = [ + ^o(Eo)dE = ^—E c - 5+1 electrons s- 1 (3.1) 
Je c o — l 

Pnth = k e J b °°E • J? (£o) dE = y^E c ~ s+2 erg s" 1 (3.2) 

The last expression in each equation is the result for a power-law electron flux distribution 
of the form &q(Eq) = A ■ Eq S . The constant Kg = 1.60 x 10~ 9 is the conversion from keV 
to erg. E c is a low-energy cutoff to the electron flux distribution. These expressions are 
valid and finite for 8 > 2 and E c > 0. We call this form of low-energy cutoff a sharp low- 
energy cutoff. An electron distribution that continues below a transition energy E c that has a 
positive slope, is flat, or in general has a spectral index 8i ow < 1 also provides finite electron 
and energy fluxes, but these fluxes are somewhat higher than those associated with the sharp 
low-energy cutoff. 

For this single-power-law electron flux distribution with a sharp low-energy cutoff, the 
non-thermal power (erg s ), and ultimately the non-thermal energy (erg), from the power- 
law electron flux distribution dep ends on only t hree parameters: 8, A, and E c . Observations 
indicate that 8 is gr eater than 2 foennislll985l : iLin & Schwartzll 19871 : Iwinglee et alll99ll : 
iHolman et alJ2 003). Hence, were E c = 0, the integral would yield an infinite value, a decid- 
edly unphysical result! Therefore, the power-law electron distribution cannot extend all the 
way to zero energy with the same or steeper slope, and some form of low-energy cutoff in 
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the accelerated electron spectrum must be present. As we will see, the determination of the 
energy at which this cutoff occurs is not a straightforward process, but it is the single most 
important parameter to determine (as the other two are generally more straightforward to 
determine - see Section[2]and Kontar et al. 2011). For example, with 5=4 (typical during 
the peak time of strong flares), a factor of 2 error in E c yields a factor of 4 error in />„,/,. For 
larger 8 (as found in small flares, or rise/decay phases of large flares), such an error quickly 
leads to an order of magnitude (or even greater) difference in the injected power P n/ /, and in 
the total energy in the non-thermal electrons accelerated during the flare! 



3.2 Why is the low-energy cutoff difficult to determine? 




1 10 100 

Photon energy e [keV] 



Fig. 3.1 Typical full-Sun flare spectrum. Dashed: Nonthermal thick-target spectrum from an accelerated 
electron distribution with 5=4, and a low-energy cutoff of 20 keV. Dotted: Thermal spectrum, from a plasma 
with temperature T = 20 MK and emission measure EM = 10 49 cm~ 3 . Solid: Total radiated spectrum. The 
multiple peaks in the thermal spectrum are from spectral lines, as observed by an instrument with ~ 1 keV 
spectral resolution. 

The essence of the problem in many flare spectra is summarized in Figure [3D the non- 
thermal power-law is well-observed above ~20 keV, but any revealing features that it might 
possess at lower energies, such as a low-energy cutoff, are masked by the thermal emission. 

Even if a spectrum does show a flattening at low energies that could be the result of 
a low-energy cutoff, other mechanisms that could produce the flattening must be ruled out 
(see Section 13.4b . The low-energy cutoff has the characteristic feature, determined by the 
photon energy dependence of the bremsstrahlung cross-section (see Equation 12. 9b, that the 
X-ray spectrum eventually approaches a spectral index of y w 1 at low energies (cf. lHolmanl 
120031) . It is currently impossible, however, to observe a flare spectrum to low enough photon 
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energies to see that it does indeed become this flat. Generally we can only hope to rale out 
the other mechanisms based on additional data and detailed spectral fits. 



3.3 What is the shape of the low-energy cutoff, and how does it impact the photon 
spectrum and P„//,? 

Bremsstrahlung photon spectra are obtained from convolution integrals over the electron 
flux distribution (equations 12.11 and 12.8b . Hence, featu res in an electron distribution are 
smoothed out in the resulting photon spectrum (see, e.g jBrown et al.|[2006h . 




Fig. 3.2 Different shapes of low-energy cutoff in the injected electron distribution (left) lead to slightly 
different photon spectra (right). The cutoff/turnover electron energy is E c =20 keV. The thin c urve in the righ t 
panel demonstrates how the cutoff can be masked by emission from thermal plasma. See also Holman ( 2003) 
for a thorough discussion of bremsstrahlung spectra generated from electron power-laws with cutoff. 



As can be seen in Figure [3~2l both a sharp cutoff at E c and a "turnover" (defined here to 
be a constant Fq(E) below E c , a "plateau") in the injected electron distribution lead to similar 
thick-target photon spectra. This subtle difference is difficult to discriminate observationally, 
and the problem is compounded by the dominance of the thermal component at low energies. 

A sharp cutoff would lead to plasma instabilities that should theoretically flatten the 
distribution around and below the cutoff within microseconds (see Section[6]l. On the other 
hand, the electron flux distribution below the cutoff must be flatter than E , as demon- 
strated by Eauation l3.il or the total electron number flux would be infinite. Having a constant 
value for the distribution below E c (turnover case) seems like a reason able middle ground 
and a pproximates a quasilinearly relaxed electron distribution (Sectionlrjl lKrall & Trivelpiecd 
1 19731 . Chapter 10). Coulomb collisional losses, on the other hand, yield an electron distri- 
bution that increases linearly at low energies (see Figure [3731 . leading to a photon spectrum 
between the sharp cutoff case and the turnover case. 
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Notice that the photon spectra actually flatten gradually to the spectral index of 1 at low 
energies from the spectral index of y = 8 + 1 at E c and higher energies. Below E c , it is not 
a power-law. Fitting a double power-law model photon spectrum, and using the break (i.e., 
kink) energy as the low-energy c utoff typically leads to a large error in E c (e.g.. lGan et al.l 
l200ll : ISaint-Hilaire & Benzll2005h. and hence to an even la rger error in P,,//,. 

In terms of the energetics. Isaint-Hilaire & Bend d2005h have shown that the choice of 
an exact shape for the low-energy cutoff as a model is not dramatically important. For a 
fixed cutoff energy, from Equation 13.21 it can be shown that the ratio of the power in the 
turnover model to the power in the sharp cutoff model without the flat component below 
the cutoff energy is 8 /2. In obtaining spectral fits, however, th e turnover model gives higher 
cutoff energies than the sharp cutoff model. Using simulations. ISaint-Hilaire & Bena (12005) 
found that assuming either a sharp cutoff model or a turnover model led to differences in 
Pnth typically less than ~20%. Hence, the sharp cutoff, being the simplest, is the model of 
choice for computing flare energetics. Nevertheless, knowing the shape of the low-energy 
cutoff would not only yield more accurate non-thermal energy estimates, but would be a 
source of information on the acceleration mechanism and/or propagation effects. 




Electron energy [keV] Electron energy [keV] 

Fig. 3.3 The four plots show the Coulomb-collisional evolution with column density of an injected electron 
distribution (thick, solid line). For the simple power-law case (upper left), the low-energy end of the distri- 
bution becomes linear, and the peak of the distribution is found at E pea it = E t /\/S, where S is the injected 
distribution power-law spectral index (5=4 in the plots), and E t = \/2K-N, is the initial energy that electrons 
must possess in order not to be fully stopped by a column density N t (Equation 12.41 . When a low-energy 
cutoff is pres ent, the peak of the distribution is seen to first decrease in energy until E t exceeds the cutoff 
energy (from Saint-Hilaire 2005). 



Spectral inversion methods have recently been devel oped for deducing the mean electron 
flux d istribution (Equation I2.5l l from X-ray spectra djohns & Linll 19921 : iBrown et al.ll2003L 
2006). A spectral "dip" has been found just above the presumed thermal component in some 
deduced mean electron flux distributions that may be associated with a low-energy cutoff 
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(e.g.. IPiana et al.ll2003h . In the collisional thick-target model, the slope of the high-energy 
"wall" of this dip should be linear or flatter, with a linear slope indicating the absence of 
emitting electron s in the injected electron distribution at the energies displaying this slope. 
iKontar & Browr] |2006a) have found evidence for slopes that are steeper than linear, but 
their spectra were not corrected for photospheric albedo (see Section [3~5l . Finding and 
understanding these dips is a crucial element for gaining an und erstanding of the low-energy 
prop erties of flare electron distributions (see lKontar et aljfeoi ll) . 

lEmslid l l2003l) has pointed out that the non-thermal electron distribution could seam- 
lessly merge i nto the thermal distr ibution, removing the need for a low-energy cutoff. As 
was shown bv lHolman etail ( t2003h for SOL2002-07-23T00:35 (X4.8), however, merger of 
the electron distribution into the typically derived ~ 10-30 MK thermal flare plasma gener- 
ally implies an exceptionally high energy in non-thermal electrons. Thus, for a more likely 
energy content, a higher low-energy cutoff or a hotter plasma would need to be present in 
the target region. Any emission from this additional "hot core," because of its much lower 
emission measure, is likely to be masked by the usual ~ 10-30 MK thermal emission. This 
merger of the non-thermal electron distribution into the thermal tail in the target region does 
not remove the need for a low-energy cutoff in the electrons that escape the acceleration 
region, however. 

This section has dealt with the shape of the low-energy cutoff under the assumptions that 
the X-ray photon spectra are not altered by other mechanisms and that the bremsstrahlung 
emission is isotropic. The next section lists the important caveats to these assumptions, and 
their possible influence in the determination of the low-energy cutoff to the electron flux 
distribution. 

3.4 Important caveats 

As previously discussed, apparently minor features in the bremsstrahlung photon spectrum 
can have substantial implications for the mean electron flux and, consequently, the injected 
electron distribution. This means that unknown or poorly-understood processes that alter 
the injected electron distribution (propagation effects, for example) or the photon spectrum 
(including instrumental effects) can lead to significant errors in the determination of the 
low-energy cutoff. Known processes that affect the determination of the low-energy cutoff 
are enumerated below. 

1. Detector pulse pileup effects dSmith et al.ll2002h . if not properly corrected for, can in- 
troduce a flattening of the spectrum toward lower energies that simulates the flattening 
resulting from a low-energy cutoff. 

2. The contribution of Compton back-scattered photons (photospheric albedo) to the mea- 
sured X-ray spectrum can sim ulate the spectral flattening produced by a low-energy 
cutoff. IKasparova et al.l ((2005) have shown that the dip in a spectrum from SOL2002- 
08-20T08:25 (M3.4) becomes st atistically insignificant when the spectrum i s corrected 
for photospheric albedo (also see lKontar et al.ll201 lh . IKasparova et al.l [2007) show that 
spectra in the 15-20 keV energy band tend to be flatter near disk center when albedo 
from isotropically emitted photons is not taken into account, further demonstrating the 
importance of correcting for photospheric albedo. 

3. The assumed differential cross-se ction and electron energy lo ss rate can influence the 
results (for a discussion of this, see Saint-Hilaire & Benz 2005). In some circumstances, 
a contribution from recombination radiation may significantly change the results (Brown 
et al. 2010; also see Kontar et al., 201 1). 
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4. Anisotropics in the electron beam directivity and the bremsstrahlung d ifferential cross- 
section can significantly alter the X-ray spectrum dMassone et al.ll2004T) . 

5. Non-uniform target ionization (the fact that the chromosphere's ionization state varies 
with depth, see Section[4]l can introduce a spectral break that may be confused with the 
break associated with a low-energy cutoff. 

6. Energy losses associated with a return current produce a low-energy flattening of the 
X-ray spectrum (Section [5}. This is a low-energy "cutoff" in the electron distribution 
injected into the thick target, but it is produced between the acceleration region and the 
emitting source region. 

7. A non-power-law distribution of injected electrons or significant evolution of the in- 
jected electron distribution during the observational integration time could affect the 
deduced value of the low-energy cutoff. 

For all the above reasons, the value of the low-energy cutoff in the injected electron flux 
distribution has not been determined with any degree of certainty except perhaps in a few 
special cases. Even less is known about the shape of the low-energy cutoff. The consensus 
in the solar physics community for now is to assume the simplest case, a sharp low-energy 
cutoff. Existing studies, presented in the next section, tend to support the adequacy of this 
assumption for the purposes of estimating the total power and energy in the accelerated 
electrons. 



3.5 Determinations of E c and electron energy content from flare data 

Before RHESSI, instruments did not cover well (if at all) the ~ 10^10 keV photon ener- 
gies where the transition from thermal emission to non-thermal emission usually occurs. 
Researchers typically assumed an arbitrary low-energy cutoff at a value at or below the in- 
strument's observing range (one would talk of the "injected pow er in electrons ab ove E c 
keV" instead of the total non-thermal power P„r/,). An exception is lNitta et al.l Jl990h . They 
argued that spectral flattening observed in two flares wi th the Sola r Maxi mum Mission and 
Hinotori indicated a cutoff energy of >50 keV. Also, iGan et alJ foOOlh interpreted spec- 
tral breaks at ~80 keV in Compton Gamma-Ray Observatory (CGRO) flare spectra as the 
low-energy cutoff in estimating flare energetics, resulting in rather small values for the non- 
thermal energy in the analyzed flares. The relatively low resolution of the spectra from these 
instruments prevented the quantitative evaluation of any spectral flattening toward lower 
energies, however. 

The only high-resolutio n flare spectral data before the launch of RHESSI was the bal- 
loon data of Lin et al.l d 198 lh for SOL198 0-06-27T16:17 (M6.7) alo ng with -25 microflares 
observed during the same balloon flight. Ben ka & Holmanl dl994l) applied a direct electric 
field electron acceleration model to the SOL1980-06-27 flare data. They derived, along with 
other model-related parameters, the time evolution of the critical energy above which run- 
away acceleration occurs - the model equivalent to the low-energy cutoff. They found this 
critical energy to range from ~ 20^10 keV. 

It is now possible in most cases to obtain a meaningful upper limit on E c , thanks to 
RHESSI's high-spec t ral-resolution cover age o f the 10—40 keV energy ran ge and beyond. 
IHolman et al]j2003h . lEmslie et al.l ll2004l) . and lSaint-Hilaire & Benzl [2005b . in determining 
the low-energy cutoff, obtained the "highest value for E c that still fits the data." In many solar 
flare spectra, because of the dominance of radiation from thermal plasma at low energies, 
a range of values for E c fit the data equally well, up to a certain critical energy, above 
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which the x 2 goodness-of-fit parameter becomes unacceptably large. The low-energy cutoff 
is taken to be equal to this critical value. This upper limit on the cutoff energy results in a 
lower limit for the non-thermal power and energy. The results obtained for the maximum 
value of E c were typically in the 15 — 45 keV range, although late in the development of 
SOL2002-07-23T00:35 (X4.8) some values as high as ~80 keV were obtained for E c . The 
minimum non-thermal energies thus determined were comparable to or somewhat larger 
than the calculated thermal energies. 
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Fig. 3.4 RHESSI spatially integrated spectra in four time intervals during SOL2002-04-15T03:55 (Ml. 2). 
(a) Spectrum at 23:06:20-23:06:40 UT (early rise phase), (b) Spectrum at 23:09:00-23:09:20 UT (just before 
impulsive phase), (c) Spectrum at 23:10:00-23:10:20 UT (soon after the impulsive rise), (d ) Spectrum at 
23:11:00-23:11:20 UT (at the hard X-ray peak). The plus signs with error bars represent the spectral data. 
The lines represent model spectral fits: the dashed lines are non-thermal thick-target bremsstra hlung, the 
dotted lines are thermal bremsstrahlung, and the solid lines are the summation of the two (from [Sui et al l 
120053) . 



On e of the best determinations of the low-energy cutoff so far was obtained bv lSui et al] 
d2005d) . They complemented the spatially-integrated spectral data for the SOL2002-04- 
15T03:55 (Ml. 2) limb flare with imaging and lightcurve information. Four spectra from 
this flare are shown in Figure 13.41 The earliest spectrum, before the impulsive rise of the 
higher energy X-rays, was well fitted with an isothermal model. The last spectrum, from 
the time of the hard X-ray peak, clearly shows a thermal component below ~20 keV. Of 
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particular interest is the second spectrum, showing both thermal and non-thermal fit com- 
ponents. As a consequence of the flattening of the isothermal component at low energies, 
the low-energy cutoff to the non-thermal component cannot extend to arbitrarily low ener- 
gies without exceeding the observed emission. This places a tight constraint on the value 
of the low-energy cutoff. The additional requirement that the time evolution of the derived 
temperature and emission measure of the thermal component be smooth and continuous 
throughout the flare constrains the value at other times. Applying the collisional thick-target 
model with a power-law distribution of injected electrons, they found the best cutoff value to 
be E c = 24±2 keV (roughly constant throughout the flare). The energy associated with these 
non-thermal electrons was found to be comparable to the peak energy in the X-ray-emitting 
thermal plasma, but an order of m agnitude greater than the kinetic energy of the associ- 
ated coronal mass ejection (CME) dSui et al]|2005bl) . This contrasts with results obtained 
for large flares, where the minimum energy in non-therm al electrons is typic ally found to be 
less than or on the order of the energy in the CME (e.g., lEmslie et al. 2004). 

The importance of correcting for the distortion of spectra by albedo was revealed by a 
search for low-energy cu toffs in a sample of 177 flares with relatively flat spectra (y < 4) 
between 15 and 20 keV ( Kont ar et al . 2008a). Spectra can be significantly flattened by the 
presence of albedo photons in this energy range. The X-ray spectra, integrated over the 
duration of the impulsive phase of the flares, were inverted to obtain the corresponding 
mean electron flux distributions. Eighteen of the flares showed significant dips in the mean 
electron flux distribution in the 13-19 keV electron energy range that might be associated 
with a low-energy cutoff (see Section [3~3l l. However, when the X-ray spectra were corrected 
for albedo from isotropically emitted X-rays, all of the dips disappeared. Therefore, the 
authors concluded that none of these flare electron distributions had a low-energy cutoff 
above 12 keV, the lowest electron energy in their analysis. 

Low-energy cutoffs were i dentified in the spectra of a sample of early impulsive flares 
observed by RHESSI in 2002 dSui et al.l l2007r l . Early impulsive flares are flares in which 
the >25 keV hard X-ray flux increase is delayed by less than 30 s after the flux increase at 
lower energies. The pre-impulsive-phase heating of plasma to X-ray-emitting temperatures 
is minimal in these flares, allowing the nonthermal part of the spectrum to be observed 
to lower energies. In the sample of 33 flares, 9 showed spectral flattening at low energies 
in spectra obtained throughout the duration of each flare with a 4 s integration time. After 
correcting for the albedo from isotropically emitted X-rays, the flattening in 3 of the 9 flares, 
all near Sun center, disappeared. The flattening that persisted in the remaining 6 flares was 
consistent with that produced by a low-energy cutoff. The values derived for the low-energy 
cutoff ranged from 15 to 50 keV. The authors found the evolution of the spectral break and 
the corresponding low-energy cutoff in these flares to be correlated with the non-thermal 
hard X-ray flux. Further studies are needed to assess the significance of this correlation. 

Low-energy cutoffs with values excee ding 100 keV were id entified in the spectra of the 
large flare SOL2005-01-19T08:22 (X1.3) dWarmuth et alj2009l) . The hard X-ray light curve 
of this flare consisted of multiple peaks that have been interpreted a s quasi-periodic oscilla - 
tions driven by either magnetoacoustic oscillations in a nearby loop (Nakariakov et al . 2006) 
or by super-Alfvenic beams in the vicinity of the reconnection region dOfrnan & Sul2 006). 
The high low-energy cutoffs were found in the last major peak of the series of hard X-ray 
peaks. Unlike the earlier peaks, this peak was also unusual in that it was not accompanied 
by the Neupert effect (see Section [H3J, consistent with the high values of the low-energy 
cutoff, and it exhibited soft-hard-harder rather than soft-hard-soft spectral evolution (see 
Section |9~T| >. A change in the character of the observed radio emission and movement of one 
of the two hard X-ray footpoints into a region of stronger photospheric magnetic field were 
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also observed at the time of this peak. These changes suggest a strong connection between 
large-scale flare evolution and electron acceleration. 



4 Nonuniform ionization in the thick-target region 

In the interpretation of hard X-ray (HXR) spectra in terms of the thick-target model, one 
effect which has been largely ignored until recently is that of vary ing ionization along the 
path of the thick-target beam. As first discussed by IBrownl dl973l) . the decrease of ioniza- 
tion with depth in the solar atmosphere reduces long-range collisional energy losses. This 
enhances the HXR bremsstrahlung efficiency there, elevating the high energy end of the 
HXR spectrum by a factor of up to 2.8 above that for a fully ionized target. The net re- 
sult is that a power-law electron spectrum of index 8 produces a photon spectrum of index 
7 = 8 — 1 at low and high energies (see Equation 12. 12b . but with y < 8 — 1 in between. The 
upward break, where the spectrum begins to flatten toward higher energies, occurs at fairly 
low energies, probably masked in measured spectra by the tail of the thermal component. 
The downward knee, where the spectrum steepens again to y = 8 — 1, occurs in the few 
deka-keV range, depending on the column depth of the transition zone. Thus, the measured 
X-ray spectrum may show a flattening similar to that expected for a low-energy cutoff in the 
electron distribution. 



4.1 Electron energy losses and X-ray emission in a nonuniformly ionized plasma 

The collisional energy-loss cross-section Q C {E) is dependent on the ionization of the back- 
ground medium. Flare-accelerated electron beams can propagate in the fully io nized corona 
as wel l as i n the par t ially io nized transition region and chromosphere.. Following Havakawa & Kitao 

and IBrownl ( 119730 . the cross-section Q C {E) can be written for a hydrogen plasma 
ionization fraction x 

271 e A 2tte' 
Q C {E) = -g- (xA ee + ( 1 - x)A eH ) = -p-A (x + X), (4. 1) 

where e is the electronic charge, A ee the electron-electron logarithm for fully ionized media 
and A e H is an effective Coulomb logarithm for electron-hydrogen atom collisions. Numeri- 
cally A ee ~ 20 and A f # ~ 7. 1, so A = A ee — A e n — 12.9 and X = A e n /A ~ 0.55. 

Then, in a hydrogen target of ionization level x(N) at column density N(z), the energy 
loss equation for electron energy E is (cf. Equation 12.2b 

rf£ = _27C^A {x+xm = _? il+x(m (42) 

where K' = 2ne 4 A = (A /A ee )K ~ 0.65£. 

The energy loss of a given particle with initial energy Eq depends on the column density 
N(z) = §Qn(z!)dz' , so the electron energy at a given distance z from the injection site can be 
written E 2 = Eq - 2K'M(N(z)) (cf. EquationHS, where 



N(z) 

M{N{z))= J (X+x(N'))dN' 
o 



(4.3) 
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Fig. 4.1 Photon spectrum residuals, normalized by the statistical error for the spectral fit, for the time interval 
00:30:00 - 00:30:20 UT, 2002-July-23, for {upper panel) an isothermal Maxwellian plus a power-law and 
(lower panel) a n isothermal Ma xwellian plus the nonuniform ionization spectrum with 5 = 4.24 and E, = 
53 keV (from Kontar et al. 2003). 



is the "effective" ionization-weighted collisional column density. 

The fractional atmospheric ionization x as a function of column density N (cm -2 ) changes 
from 1 to near over a small spatial range in the solar atmosphere. Therefore, to lowest 
order, x{N) can be approximated by a step function x(N) = 1 for N < N*, and x(N) = 
for N > N*. This gives M(N) = (A + l)N for N < TV* and M(N) =N*+ AiV fo r N > N*. 
Electrons injected into the target with energies less than E„ = \J2K'(X + 1)2V* = \/2KN* 
experience energy losses and emit X-rays in the fully ionized plasma with x = 1, as in the 
standard thick-target model. Electrons injected with energies higher than E. t lose part of 
their energy and partially emit X-rays in the un-ionized (x = 0), or, more generally, partially 
ionized plasma. 

We can deduce the properties of the X-ray spectrum by substituting Equation 14.21 into 
Equation l2.7l (with clN = nvdt) and comparing Ithick(^) from Equation 12 . 6 1 with 7, ;„>i- ( e ) from 
Equation 12.81 We see that for the nonuniformly ionized case the denominator in the inner 
integral now contains A +x(N) and K is replaced with K' . In the step-function model for 
x{N), photon energies greater than or equal to = are emitted by electrons in the un- 
ionized plasma with E > Est ■ Since A +x(N) has the constant value A , the thick-target power- 
law spectrum is obtained (for injected power-law spectrum), but the numerical coefficient 
contains K'X = 27ie 4 A e n instead of K. At photon energies far enough below that the 
contribution from electrons with E > E t is negligible, A +*(iV) = A + 1 and the numerical 
coefficient contains (A + \)K' = K. The usual thick-target spectral shape and numerical 
coefficient are recovered. The ratio of the amplitude of the high-energy power-law spectrum 
to the low-energy power-law spectrum is (A + 1)/A ~ 2.8. The photon energy £*(keV) ~ 
2.3 x 10~ 9 ^/N* (cm -2 ), between where the photon spectrum flattens below the high-energy 
power law and above the low-energy power law, determines the value of the column density 
where the plasma ionization fraction drops from 1 to 0. 



4.2 Application to flare X-ray spectra 



The step-function nonuniform ionization model was used bv iKontar et ail d2002l . l2003h to 
fit photon spectra from five flares. They assume a single power-law distribution of injected 



20 



Holman et al. 



electrons with power-law index 8 and approximate the bremsstrahlung cross-section with 
the Kramers cross-section. First, they fit the spectra to the sum of a thermal Maxwellian 
at a sin gle temperature T plus a single power law of index 7. For SOL2002-07-23T00:35 
(X4.8) kontar et al.ll2003h they limit themselves to deviations from a power law in the non- 
thermal component of the spectrum above ~40 keV. The top panel of Figure 14.11 shows an 
example of such deviations, which represent significant deviations from the power-law fit. 
These deviations are much reduced by replacing the power law with the spectrum from the 
nonuniform ionization model, with the minimum rms residuals obtained for values of 8 = 
4.24 and = 53 keV (Figure |4~T1 bottom panel). The corresponding minimum (reduced) 
X 2 value obtained for the best fit to the X-ray spectrum (10-130 keV) dropped from 1.4 for 
the power-law fit to 0.8 for the nonuniform ionization fit. There are still significant residuals 
present in the range from 10 to 30 keV; these might be due to photospheric albedo or the 
assumption of a single-temperature thermal component. 

By assuming that the main spectral feature observed in a hard X-ray spectrum is due 
to the increased bremsstrahlung efficiency of the un-ionized chromosphere, allowance for 
nonuniform target ionization offers an elegant direct explanation for the shape of the ob- 
served hard X-ray spectrum and provides a measure of the location of the transition re- 
gion. Table 14. II sho ws the best fit parameters derived for the four flare spectra analyzed by 
iKontar et alJ (120021) . The last column shows the ratio of the minimum x 2 value obtained 
from the nonuniform ionization fit to the minimum x 2 value obtained from a uniform ion- 
ization (single power-law) fit to the non-isothermal part of the spectrum. The nonuniform 
ionization model fits clearly provide substantially better fits than single power-law fits. 



Table 4.1 Best fit nonuniformly ionized target model parameters for a single power-law .j Pq(Eq), the equiv- 
alent N, (energy range 20-100 keV), and the ratio of xlonunil %lni (from lKontar et all2002ri 



Date 


Time, UT 


ytr(keV) 


S 


E, (keV) 


N t 


(cm 2 ) 


7 2 -It 2 ■ 


20 Feb 2002 


11:06 


1.47 


5.29 


37.4 


2.7 


<10 2u 


0.032 


17 Mar 2002 


19:27 


1.27 


4.99 


24.4 


1.1 xlO 20 


0.047 


31 May 2002 


00:06 


2.02 


4.15 


56.2 


6.1 


<10 20 


0.041 


1 Jun 2002 


03:53 


1.45 


4.46 


21.0 


8.4 


<10 19 


0.055 



Values of the fit parameters kT (keV), 8 and E* as a function of time for SOL2002-07- 
23T00:35 (X4.8) together with the c orresponding value of AT* (cm -2 ) ~ 1.9 X 10 17 £*(keV) 2 
were obtained bv lKontar et alJ I I2003I) . The results (Figure 14721 demonstrate that the thermal 
plasma temperature rises quickly to a value ~ 3 keV and decreases fairly slowly thereafter. 
The injected electron flux spectral index 8 follows a general "soft-hard-soft" t rend and qual- 
itativ ely agrees with the time history of the simple best-fit power-law index 7 ( iHolman et al.l 
120031) . £* rises quickly during the first minute or so from ~40 keV to ~70 keV near the 
flare peak and thereafter declines rather slowly. The corresponding values of N* are ~2- 
5 xlO 20 cm" 2 . 

The essential results of these studies are that (1) for a single power-law electron injection 
spectrum, the expression for bremsstrahlung emission from a nonuniformly-ionized target 
provides a significantly better fit to observed spectra than the expression for a uniform target; 
and (2) the value of E* (and correspondingly 2V*) varies with time. 

An upper limit on the d egree o f spect ral flattening A 7 that can result from nonuniform 
ionization was derived by ISu et all d2009h . They applied this upper limit to spectra from a 
sample of 20 flares observed by RHESSI in the period 2002 through 2004. They found that 
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Fig. 4.2 Variation of kT , 8, £», and N* throughout SOL2002-07-23T00:35 O C4.8) iKontar et aljfc tXB). 
The variation of o ther parameters, such as emission measure, can be found in iHolman^HTl . 1 2003) and 
ICaspi & LinU201Cl) . 



15 of the 20 flare spectra required a downward spectral break at low energies and for each of 
these 15 spectra derived the difference Ay of the best- fit power-law spectral indices above 
and below the break. A Monte Carlo method was used to determine the 95% confidence 
interval for each of the derived values of Ay. Taking the value of A y to be incompatible with 
nonuniform ion ization if the 95% confidence interval fell above the derived upper limit, 
ISu et alJ (2009) found that six of the flare spectra could not be explained by nonuniform 
ionization alone. Thus, for these six flares some other cause such as a low-energy cutoff or 
return-current-associated energy losses (Section|5]l must be at least partially responsible for 
the spectral flattening. 



5 Return current losses 

The thick-target model assumes that a beam of electrons is injected at the top of a loop and 
"precipitates" downwards in the solar atmosphere. Unless accompanied by an equal flux of 
positively charged particles, these electrons constitute a current and must create a significant 



22 



Holman et al. 



self-induced electric field that in turn drives a co-spatial return current for compens ation 
dHovng et al.ll976l : lKnight & SturroclJl1^ : lEmsliell980l : lDTakonov & Somovll988h . The 
return current consists of ambient electrons, plus any primary electrons that have scattered 
back into the upward direction. By this means we have a full electric circuit of precipitating 
and returning electrons that keeps the whole system neutral and the electron beam stable 
against being pinched off by the self-generated magnetic field required by Ampere's law for 
an unneutralized beam current. However, the self-induced electric field results in a potential 
drop along the path of the electron beam that decelerates and, therefore, removes energy 
from the beam electrons. 

5.1 The return current electric field 

The in itial formation of the beam/retum-current system has been studied by Ivan den Oordl 
dl990h and references therein. We assume here that the system has time to reach a quasi- 
steady state. Van den Oord finds this time scale to be on the order of the thermal electron-ion 
collision time. This time scale is typically less than or much less than one second, depend- 
ing on the temperature and density of the ambient plasma. In numerical simulations by 
I Siverskv & Zharkoval d2009h . times to reach a steady state after injection ranged from 0.07 s 
to 0.2 s, depending on the initial beam parameters. 

The self-induced electric field strength at a given location z along the beam and the flare 
loop, £(z), is determined by the current density associated with the electron beam, j(z), 
and the local conductivity of the loop plasma, cr(z), through Ohm's law: £{z) = j(z)/o(z). 
Relating the current density to the density distribution function of the precipitating electrons, 
f(z,E, 9), where E is the electron energy and 6 is the electron pitch angle, gives 

S(z) = f j ' f( z ,E,B)y/EndEdn. (5.1) 

0(Z) Jm e J J 


Here /i is the cosine of the pitch angle and e and m e are the electron charge and mass, 
respectively. The self-induced electric field strength S(z) depends on the local distribution 
of the beam electrons, which in turn depends on the electric field already experienced by 
the beam as well as any Coulomb energy losses and pitch-angle scattering that may have 
significantly altered the beam. It also depends on the local plasma temperature (and, to a 
lesser extent, density) through c(z), which can, in turn, be altered by the interaction of the 
beam with the loop plasma (i.e., local heating and "chromospheric evaporation"). Therefore, 
determination of the self-induced electric field and its impact on the precipitating electrons 
generally requires self-consistent modelin g of the coupled beam/plasma system. 

Such models have been computed bv lZharkova & Gordovskv v (200 5L l2006h . They nu- 
merically integrate the time-dependent Fokker-Planck equation to obtain the self-induced 
electric field strength and electron distribution function along a model flare loop. The in- 
jected electron beam was assumed to have a single power-law energy distribution in the 
energy range from E[ OH , = 8 keV to E upp = 384 keV and a normal (Gaussian) distribution in 
pitch-angle cosine fi with half- width dispersion Apt = 0.2 about fi = 1, 

The model computations show that the strength of the self-induced electric field is nearly 
constant at upper coronal levels and rapidly decreases with depth (column density) in the 
lower corona and transition region. The rapidity of the decrease depends on the beam flux 
spectral index. It is steeper for softer beams (5=5-7) than for harder ones (5=3). The strength 
of the electric field is higher for a higher injected beam energy flux density (erg cm~ 2 s _1 ), 
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Photon Energy (keV) 

Fig. 5.1 (a) Photon spectra computed from full kinetic solutions including return current losses and colli- 
sional losses and scattering. The top spectrum is for an injected single-power-law electron flux distribution 
between 8 keV and 384 keV with an index of 5 = 3, and the bottom spectrum is for 5 = 7. The injected 
electron energy flux density is 10 s erg cm~ 2 s -1 . (b) Same as (a), but for an injected energy flux density of 
10 12 erg cm -2 s _I . The tangent lines at 20 and 100 keV demonstrate the determination of the low-energy and 
high-energy power-law spectral indices y/ OK . and y/„ s /, . (c) The photon spectral indices yj ow (dashed lines) and 
Thigh (solid lines) vs. 5 for an injected energy flux density of 10 8 (squares), 10 10 (circles), and 10 12 erg cm -2 
s (crosses), (d) %i s h — "flaw vs. the log of the injected electron energy flu x density for S equal to 3 (bottom 
curve, squares), 5 (middle curve, circles), and 7 (top curve, triangles) (from Zharkova & Gordovskyy 2006). 

and the distance from the injection point over which the electric field strength is highest (and 
nearly constant) decreases with increasing beam flux density. 



5.2 Impact on hard X-ray spectra 

Deceleration of the precipitating beam by the electric field most significantly affects the 
lower energy electrons (<100 keV), since the fraction of the original particle energy lost 
to the electric field is greater for lower energy electrons. This leads to flattening of the 
electron distribution function towards the lower energies and, therefore, flattening of the 
photon spectrum. 

Photon spectra computed from kinetic solutions that include return current energy losses 
and collisional energy losses and scattering are shown in Figure IBTTl (a) and (b). Low- and 
high-energy spectral indices and their dependence on the power-law index of the injected 
electron distribution and on the injected beam energy flux density are shown in Figure |5TI 
(c) and (d). The difference between the high-energy and low-energy spectral indices is seen 
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to increase with both the beam energy flux density and the injected electron power-law index 
8. The low-energy index is found to be less than 2 for 8 as high as 5 when the energy flux 
density is as high as 10 12 erg cm~ 2 s -1 . 



5.3 Observational evidence for the presence of the return current 

We have seen that return current energy losses can introduce curvature into a spectrum, 
possibly explaining the "break" often seen in observed flare X-ray spectra. A difficulty in 
directly testing this explanation is that the thick-target model provides the power (energy 
flux) in the electron beam (erg s -1 ), but not the energy flux density (erg cm~ 2 s _1 ). X-ray 
images provide information about the area of the target, but this is typically an upper limit 
on the area. Even if the source area does appear to be well determined, the electron beam 
can be filamented so that it does not fill the entire area (the filling factor is less than 1). 
Also, if only an upper limit on the low-energy cutoff to the electron distribution is known, as 
described in Section [331 the energy flux density may be higher. Therefore, the observations 
typically only give a lower limit on the beam energy flux density. 

The non-thermal hard X-ray flux is proportional to the electron beam flux density, but 
the retu rn-current ener gy losses are also proportional to the beam flux density. As a conse- 
quence, lEmsliel dl980h concluded that the flux density of the non-ther mal X-ray emission 
from a flare cannot exceed on the order of 10~ 15 cm~ 2 s _1 above 20 ke V. I Alexander & Daoul 
J2007h have deduced the photon flux density from non-thermal electrons in a sample of 10 
flares ranging from GOES class Ml. 8 to X17. They find that the non-thermal photon flux 
density does not monotonically increase with the thermal energy flux, but levels off (satu- 
rates) as the thermal energy flux becomes high. They argue that this saturation most likely 
results from the growing importance of return current energy losses as the electron beam flux 
increases to high values in the larger flares. They find that the highest non-thermal photon 
flux densities agree with an upper limit computed by Emslie. 

A correlation between the X-ray flux and spectral break energy was found bv lSui et"aD 

d2007l) in their study of X-ray spectra in early impulsive flares (see Section [3~5l >. They point 
out that the increasing impact of return current energy losses on higher energy electrons as 
the electron beam energy flux density increases could be an explanation for this correlation. 

iBattaelia & Bend d2008l) studied two flares with non-thermal coronal hard X-ray sources 
for which the difference between the measured photon spectral index at the footpoints and 
the spectral index of the coronal source was greater than two, the value expected for coronal 
thin-target emission and footpoint thick-target emission from a single power-law electron 
distribution (see Section [l0.2b . They argue that return-current losses between the coronal 
and footpoint source regions are most likely responsible for the large difference between the 
spectral indices. 

The return current can also affect the spectral line emission from flares. Evidence for 
the presence of the return current at the chromospheric level from obse rvations of the linear 
polari z ation of the hydrogen Ha and H/3 lines has been presented by Henoux & Karlickv 
J2003I) . Dzifcakova & Karlickvl (120081) have shown that the presence of a return current in 
the corona may have a distinguishable impact on the relative intensities of spectral lines 
emitted from the corona. 
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6 Beam-plasma and current instabilities 

Interaction of the accelerated electrons with plasma turbulence as they stream toward the 
thick-target emission region can modify the electron distribution. In this section we briefly 
discuss a likely source of such turbulence: that generated by the electron beam itself. If the 
beam is or becomes unstable to driving the growth of plasma waves, these waves can interact 
with the beam and modify its energy and/or pitch angle distribution until the instability 
is removed or the wave growth is stabilized. The return current associated with the beam 
(Section[5]l can also become unstable, resulting in greater energy loss from t he beam. Beam- 
plasma and return-current instabilities in solar flares have been reviewed bv lMelrosej ( fl990h 
and lBendj2002h . 

A sharp, low-energy cutoff to an electron beam or, more generally, a positive slope in 
the beam electron energy distribution is well-known to generate the growth of electrostatic 
plasma waves (Langmuir waves). The characteristic time scale for the growth of these waves 
is on the order of [(Nb/n)co pe ]~ , where co pe is the electron plasma frequency and Ni,/n the 
ratio of the density of unstable electrons in the beam to the plasma density. This is on the 
order of microseconds for a typical coronal loop plasma density and Nj,/n w 10~ 3 . This 
plasma instability is often referred to as the bump-on-tail instability. The result is that on 
a somewhat longer but comparable time scale electrons from the unstable region of the 
electron distribution lose energy to the waves until the sharp cutoff is flattened so that the 
distribution no longer drives the rapid growth of the waves. Therefore, the electron energy 
distribution below the low-energy cutoff is likely to rapidly become flat or nearly flat (suffi- 
ciently flat to stabili ze the instability) after the electrons escape the acceleration region (see 
Chapters 9 & 10 of iKrall & Trivelpiecelll973h . 

A recent simulation of the bump-on-tail instability for solar flare conditions, including 
Coulomb collisions and wave da mping as the electr ons propagate into an increasingly dense 
plasma, has been carried out bv lHannah et al.l ( 120091) . The authors compute the mean elec- 
tron flux distribution in their model flare atmosphere after injecting a power-law distribution 
with a sharp low-energy cutoff. They find a mean electron flux distribution with no dip (see 
Section 13.3b and a slightly negative slope below the original cutoff energy with a spectral 
index 8 between and 1 . 

A beam for which the mean electron velocity parallel to the magnetic field substantially 
exceeds the mean perpendicular velocity can drive the growth of waves that resonantly in- 
teract with the beam. When the electron gyrofrequency exceeds the plasma frequency, these 
waves are electrostatic and primarily scatter the electrons in pitch angle. Generally known 
as the ano malous Doppler reso nance instability, this instability tends to isotropize the beam 
electrons. iHolman et al.l ([1982) showed that under solar flare conditions this instability can 
grow and rapidly isotropize the beam electrons in less than a millisecond. They found that 
electrons at both the low- and high-energy ends of the distribution may remain unscattered, 
however, because of wave damping . This could result in up t o two breaks in the emitted X- 
ray spectrum. On the other hand, Vlah os & Rowland dl984l) have argued that non-thermal 
tails will form in the ambient plasma and stabilize the anomalous Doppler resonance insta- 
bility by suppressing the growth of the plasma waves. 

Electrons streaming into a converging magnetic field can develop a loss-cone distribu- 
tion, with a deficit of electrons at small pitch angles. Both classical Coulomb collisions 
and loss-cone instabilities can relax this distribution by scattering electrons into the loss 
cone or extracting en ergy from the comp onent of the electron velocities perpendicular to the 
magnetic field (e.g., lAschwander] [l990). One loss-cone instability, the electron-cyclotron 
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or gyrosynchrotron maser, produces cohe rent radiation observable at radio frequencies 
toolrnanetalJl98d : iMekose & Dulll982h . 

The return current associated with the streaming electrons becomes unstable to the ion- 
acoustic instability when its drift speed exceeds a value on the order of the ion sound speed. 
The excited ion sound waves enhance the plasma resistivity, increasing the electric field 
strength associated with the return current, the heating of the plasma by the current, and the 
energy loss from the electron beam. 

It has been argued that rapid plasma heating and particle acceleration in the corona 
should result in the expansion of hot plasma down the legs of flare lo ops at the ion soun d 
speed, confined behind a collisionless ion-acoustic conduction front dBrown et al.| [l979). 
Electrons with speeds greater than about three times the electron thermal speed would be 
able to stream ahead of the conduction front. This scenario has not been observationally 
verified, but the observational signature may be confused by the chromospheric evaporation 
produced by the high-energy particles streaming ahead of the conduction front. 

iRowland & Vlahosl l l 19851) argued that if the electron beam is unstable to beam plasma 
interactions, the return current will be carried by high-velocity electrons. This reduces the 
impact of collisio ns on the beam/return -current system and helps stabilize the system. In a 
recent simulation, iKarlickv et all J2008h have found that for current drift velocities exceed- 
ing the electron thermal speed, the return current is carried by both the primary (drifting 
thermal) current and an extended tail of high- velocity electrons. 

The evolution of the electron-beam/return-current system when the return-curre nt drift 
speed exceeds the electron thermal speed has also been simulated bv lLee et"aD j2008h . They 
find that double layers form in the return current, regions of enhanced electric field that 
further increase the energy losses of the electron beam. This, in turn, increases the high- 
est electron energy to which these losses significantly flatten the electron distribution and 
corresponding hard X-ray spectrum. 

The beam/return-current system has been simulated bv lKarlickvl ( l2009]) . with a focus on 
the role o f the Weibel instability. The W eibel instability tends to isotropize the electron dis- 
tribution. Karlicky & Kasparoval d2009t> have computed the thin-target X-ray emission from 
the evolved electron distributions for a model with a weak magnetic field and another model 
with a strong magnetic field (ratio of the electron gyrofrequency to the plasma frequency 
~0 and ~ 1, respectively). They demonstrate that in both cases the electron distribution is 
more isotropic and the directivity of the X-ray emission is lower than when the instability of 
the system is not taken into account, with the greatest isotropization occurring in the weak 
field limit. 

Although we expect plasma instabilities to affect the evolution of the electron beam, 
observationally identifying them is difficult. The bump-on-tail instability and return current 
losses both lead to a flat low-energy cutoff. So far we have not established the ability to 
observationally distinguish a flat low-energy cutoff from a sharp low-energy cutoff. The 
bump-on-tail instability may be distinguishable from return-current losses by its short time 
scale and, therefore, the short distance from the acceleration region over which it effectively 
removes the unstable positive slope from the electron energy distribution. The instabili- 
ties that isotropize the electron pitch-angle distribution may be responsible for evidence 
from albedo measuremen ts that flare electron distributions are isotropic or nearly isotropic 
dKontar & Browrj|2006bl) . 
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7 Height dependence and size of X-ray sources with energy and time 

7.1 Footpoint Sources 

Hard X-ray footpoint sources result from collisional bremsstrahlung radiated by precipitat- 
ing electrons, which produce most of the emission in the chromosphere according to the 
collisional thick-target model. Depending on the density structure in the legs of the coronal 
magnetic loop, mildly energetic electrons lose their energy in the lower corona or transi- 
tion region, while the more energetic electrons penetrate deeper into the chromosphere (see 
Equation l2.2b . 

The altitude of these hard X-ray footpoint sources could never be measured accurately 
before RHESSI, because of a lack of spatial and spectral resolution. With RHESSI, we can 
measure the centroid of the footpoint location with an accuracy of order an arcsecond for 
every photon energy in steps as small as 1 keV. For a flare near the limb (Figure |TT}. the 
centroid location translates directly into an altitude. 

lAschwanden et al.l l t2002h studied such a flare, SOL2002-02-20T1 1:07 (C7.5). The heights 
of the footpoint sources were fitted with a power-law function of the photon energy. This 
yielded altitudes h w 1000— 5000 km in the energy range e = 10 — 60 keV, progressively 
lower with higher energy, as expected from the thick-target model (Figure 17711 right frame). 

Since the stopping depth of the precipitating electrons is a function of column density, 
the integrated density along their path in the chromosphere (equation 12.4b . the measured 
height dependence of the hard X-ray ce ntroids can be inverted to yield a density model 
of the chromosphere ferown et al. Assuming the decrease in density with height 

had a power-law dependence and the plasma is fully ionized, the inversion of the RHESSI 
data in the example shown in Figure 17.11 yielded a chromospheric density model that has 
a significantly higher electron density in the h = 2000 — 5000 km range than the standard 
chromospheric models based on UV spectroscopy and hydrostatic equilibrium (VAL and 
FAL models). The RHESSI-based chromospheric density model was therefore found to 
be more consistent with the "spicular extended chromosphere," simila r to the results fro m 
sub-mm radio observations during solar eclipses carried out at Caltech tewelletalJI 19931) . 

Forward fitting RHESSI X-ray visibilities to an assumed circular Gaussian source shape, 
iKontar et al.l ((2008b) found for a limb flare the full width at half maximum (FWHM) size 
and centroid positions of hard X-ray sources as a function of photon energy with a claimed 
resolution of ~0".2. They show that the height variation of the chromospheric density and 
of the magnetic flux density can be found with a vertical resolution of ~ 150 km by mapping 
the 18 — 250 keV X-ray emission of energetic electrons propagating in the loop at chromo- 
spheric heights of 400— 1500 km. Assuming collisional losses in neutral hydrogen with an 
exponential decrease in density with height, their observations of SOL2004-01-06T06:29 
(M5.8) suggest that the density of the neutral gas is in good agreement with hydrostatic 
models with a scale height of around 140 ± 30 km. FWHM sizes of the X-ray sources de- 
crease with energy, suggesting the expansion (fanning out) of magnetic flux tubes in the 
chromosphere with height. The magnetic scale height B(z)(dB/dz)~ l is found to be on the 
order of 300 km and a strong horizontal magnetic field is associated with noticeable flux 
tube expansion at a hei ght of ^900 km. A subsequent analysis with an assumed elliptical 
Gaussian source shape llKontar et al.ll2010T) confirms these results and shows that the verti- 
cal extent of the X-ray source decrease with increasing X-ray energy. The authors find the 
vertical source sizes to be larger than expected from the thick-target model and suggest that 
a multi-threaded density structure in the chromosphere is required. The thick-target model 
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Fig. 7.1 The centroids of footpoint hard X-ray emission are marked for different photon energies between 10 
keV and 60 keV for SOL2002-02-20T1 1 :()7 (C7.5), which occurred near the solar west limb and was imaged 
with RHESSI (left panel). The altitud e h(e) as a function of en ergy e shows a systematic height decrease with 
increasing energy (right panel) (from Aschwanden et al. 2002). 

to which the results were compared, however, did not account for partial occultation of the 
X-ray sources by the solar limb. 

The flare SOL2002-02-20T11:07 (C7.5) has been reanalyzed bv lPrato et all J2009h us- 
ing both photon maps over a range of photon energies and mean electron flux maps deduced 
from RHESSI visibilities over a range of electron energies. Using source centroids computed 
from the maps and assuming an exponential decrease in density with height, they found the 
density scale height to be an order of magnitude larger than the expected chromospheric 
scale height on the quiet Sun, but consistent with the scale height in a non-static, flaring 
atmosphere. This is also consistent with the enhanced plasma densities found at ~1000- 
5000 km altitudes bv lAschwanden et all J2002h . 

If the results for th e 400-1500 km height range ( Kontar et al. 2008b) and for the ~ 1000- 
5000 km height range dAschwanden et al]l2002l : |Prato et al.ll2009h are typical of flare loops, 
they imply that the upper chromosphere and transition region respond with a non-hydrostatic, 
expanded atmosphere while the low chromosphere does not respond to the flare energy re- 
lease. These results could, of course, depend on the magnitude of the flare. More studies of 
this kind are clearly desirable, especially in coordination with observations of spectral lines 
from the chromosphere and transition region. 



7.2 Loop Sources and their Evolution 

As discussed above, footpoint sources are produced by bremsstrahlung emission in the thick- 
target chromosphere. The compactness of such sources results from the rapid increase of the 
density from the tenuous corona to the much denser chromosphere. This also gives rise 
to the compact height distribution of emission centroids at different energies as shown in 
Figure [7TI However, if the density distribution has a somewhat gradual variation, one would 
expect a more diffuse height distribution. Specifically, at some intermediate energies, we 
expect that HXR emission would appear in the legs of the loop, rather than the commonly 
observed looptop sources at low energies and footpoint sources at high energies. This has 
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Fig. 7.2 CLEAN images at 04:58:22-04:58:26 UT during the impulsive phase of SOL2003-1 1-13T05:01 
(Ml. 6). The background shows the image at 9-12 keV. The contour levels are 75% and 90% of the peak flux 
at 9-1 2 keV {looptop) , 70% and 90% at 12-18 keV (legs), and 50%, 60%, & 80% at 28^13 keV (footpoints) 
(from lLiu et alj2006h . 

been observed by RHESSI in SOL200 3-11-13T05:01 ( Ml. 6) dLiu et alj|20Q6h (Figure [7T2l 
and in SOL2002-11-28T04:37 (C1.0) dSui et al.ll2006bl) . 

To reveal more details of the energy-dependent structure of SOL2003-1 1-13T05:01 
(Ml. 6), Figures U3b -c show the X-ray emission profile along the flare loop at different 
energies for three time intervals in sequence. The high energy emission is dominated by the 
footpoints, but there is a decrease of the separation of the footpoints with decreasing energy 
and with time. At later times the profile becomes a single source, peaking at th e looptop. 
The g eneral trend suggests an increase of the plasma density in the loop with time dLiu et al.l 
2006), which can be produced by chromospheric evaporation and can give rise to progres- 
sively shorter stopping distances for electrons at a given energy. Such a density increase also 
smooths out to some extent the sharp density jump at the transition region. This results in 
the non-thermal bremsstrahlung HXRs at intermediate energies appearing in the legs of the 
loop, at higher altitudes than the footpoints, as shown in Figure PT2l 

Fr om the emission profiles in the non-thermal regimes of the photon spectra, iLiu et al.l 
(2006) derived the density distribution along the lo op, using the empir i cal fo rmula for non- 
thermal bremsstrahlung emission profiles given bv lLeach & Petrosianl Jl983l . Equation 11). 
Leach and Petrosian found that this formula closely approximates their numerical results 
for a steady-state, power-law injected electron distribution with a uniform pitch-angle dis- 
tribution, no return-current losses, and a loop with no magnetic field convergence. Since this 
formula is a function of the column density, one does not need to assume any model form of 
the density distribution (cf. lAschwanden et alj|2002h . Figure I7~4l shows the density profiles 
derived from the emission profiles in the three time intervals shown in Figure 1731 Between 



30 



Holman et al. 




Fig. 7.3 (a) Brightness profiles in different energy bands measured along a semi-circular path fit to the flaring 
loop for the time interval 04:58:00-04:58:24 UT of SOL2003-1 1-13T05:01. The vertical axis indicates the 
average photon energy (logarithmic scale) of the energy band for the profile. Representative energy bands 
(in units of keV) are labeled above the corresponding profiles. The filled circles mark the local maxima, and 
the vertical dotted lines are the average positions of the centroids of the looptop and footpoint sources, (b, 
c) Same as (a), but for 04:58:24-04:58:48 an d 04:58:48-04:5 9:12 UT, respectively. The error bars show the 
uncertainty of the corresponding profile (from Liu et al. 2006). 




Fig. 7.4 Averaged density profiles along one leg of the loop inferred from the HXR brightness profiles 
during the three time intervals in Figure 1731 The distance is measured along the leg extending from the 
centroid of the thermal looptop source at about 15 arcsec in Figure 1731 to the end of the fitted semi-circle at 
about 37 arcsec (from Liu et al. 2006). 



the first and second intervals, the density increases dramatically in the lower part of the loop, 
while the density near the looptop remains essentially unchanged. The density enhancement 
then shifts to the looptop from the second to the third interval. This indicates a mass flow 
from the chromosphere to the looptop, most likely caused by chromospheric evaporation. 
For papers studying chromosp heric evaporation using coord inated RHESSI HXR and EUV 
Doppler-shift observations, see lMilligan et al] d2006al bl) and lBrosius & Holmanl ( 120071) . 
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Fig. 7.5 RHESSI (solid lines) and GOES 1-8 A (dotted line) light curves are shown in the top panel. The 
RHESSI energy bands (from top to bottom) are 3-6, 6-12, 12-25, and 50-100 keV, with scaling factors of 5, 
1, 4, 3, and 0.5, respectively. The RHESSI and GOES integration times are 4 and 3 s, respectively. The bottom 
panel shows the distance between the 3-6 keV moving source centroids and their corresponding footpoint 
centroids located in the 25-50 keV image of the flare at the time of peak emission. The distances are plane- 
of-sky values with no correction for motions away from or toward the observer (from Sui et al. 2006b). 

The flare SOL2002-1 1-28T04:37 (C1.0) was an early impulsive flare, meaning that there 
was minimal pre-heating of plasma to X-ray-emitting temperatures prior to the appear- 
ance of impulsive hard X-ray emission (see Section [331 l. RHESSI observations of this flare 
showed coronal X-ray sources th at first moved downward and then upward along the legs 
of the flare loop dSui et al.ll2006bl) . The bottom panel of Figure [73] shows the motion of the 
sources observed in the 3-6 keV band. RHESSI and GOES light curves are shown in the top 
panel. The sources originated at the top of the flare loop and then moved downward along 
both legs of the loop until the time of peak emission at energies above 12 keV. Afterward the 
source in the northern leg of the loop was no longer observable, but the source in the south- 
ern leg moved back to the top of the loop. Its centroid location at the looptop was slightly 
but significantly lower than the centroid position at the beginning of the flare. Higher-energy 
sources showed a similar evolution, but they had lower centroid positions than their lower 
energy counterparts, again in agreement with the predictions of the thick-target model. 

The early downward source motion along the legs of the loop is a previously unobserved 
phenomenon. At this time we do not know if the occurrence is rare, or if it is simp l y rarely 
observed because of masking by the radiation from the thermal plasma. ISui et ID (2006b) 
argue that the motion results from the hardening of the X-ray spectrum, and possibly an in- 
crease in the low-energy cutoff, as the flare hard X-ray emission rises to its peak intensity. A 
flatter spectrum results in a higher mean energy of the electrons contributing to the radiation 
at a given X-ray energy. In a loop with a plasma density that increases significantly from the 
top to the footpoints, these higher energy electrons will propagate to a lower altitude in the 



32 



Holman et al. 



loop as the spectrum hardens. The softening of the spectrum after peak emission would also 
contribute to the upward motion of the source after the peak. However, at that time chromo- 
spheric evaporation would likely be increasing the density in the loop, as discussed above 
for SOL2003-1 1-13T05:01 (Ml. 6) , and thermal emission would be more important. All of 
these can contribute to an increase in the height of the centroid of the X-ray source. The 
downward motion may only occur in initially cool flare loops, i.e., early impulsive flares, 

because these loops are most likely to contain the density gradients that are required. 

In an attempt to differentiate between thermal and non-thermal X-ray emission. lXu et al.l 

(2008) modeled the size dependence with photon energy of coronal X-ray sources observed 
by RHESSI in ten M-class limb flares. They determined the one-sigma Gaussian width of the 
sources along the length of the flare loops by obtaining forward fits to the source visibilities. 
The integration times ranged from one to ten minutes and the source sizes were determined 
in up to eight energy bins ranging in energy from as low as 7 keV to as high as 30 keV. They 
found the source sizes to increase slowly with photon energy, on average as e 1 / 2 . The results 
were compared with several models for the variation of the source size with energy. The 
source size was expected to vary as e -1 / 2 for a thermal model with a constant loop density 
and a temperature that decreased with a Gaussian profile along the legs of the loop from 
a maximum temperature at the top of the loop. For the injection of a power-law electron 
flux distribution into a high-density loop so that the loop is a collisional thick target, the 
source size was expected to increase as e 2 . Neither of these models are consistent with the 
observed e 1 / 2 dependence. A hybrid thermal/non-thermal model and a non-thermal model 
with an extended acceleration region at the top of the loop were found to be consistent with 
the deduced scaling, however. The extended acceleration region was deduced to have a half- 
length in the range 10" - 18" and density in the range (1 — 5) x 10 11 cm~ 3 . We note that the 
extended acceleration region model implies a column density in the range 0.73 — 6.5 x 10 20 
cm" 2 along the half length and, from Equation 12.41 all electrons with energies less than 
somewhere in the range of 23 keV - 68 keV that traverse this half length will lose all of their 
energy to collisions. The acceleration process would therefore need to be efficient enough to 
overcome these losses. On the other hand, the 7-30 keV energy range is the range in which 
fits to spatially integrated X-ray spectra typically show a combination of both thermal and 
non-thermal bremsstrahlung emission. 

Studies of flare hard X-ray source positions and sizes as a function of photon energy and 
time hold great promise for determining the height structure of flare plasma and its evolution, 
as well as providing information about the magnetic structure of the flare loop. Such studies 
are currently in their early stages, in that they usually assume an over simplified power-law 
or exponential height distribution for the plasma and do not take into account the variation 
of the plasma ionization state with height. They also assume the simple, one-dimensional 
collisional thick-target model, without consideration of the pitch-angle distribution of the 
beam electrons or the possibility of additional energy losses to the beam (such as return- 
current losses). Given the potential for obtaining a better understanding of flare evolution, 
we look forward to the application of more sophisticated models to the flare hard X-ray data. 



8 Hard X-ray timing 

The analysis of energy-dependent time delays allows us to test theoretical models of physi- 
cal time scales and their scaling laws with energy. In the wavelength domain of hard X-rays, 
there are at least three physical processes known in the observatio n of solar flares that lead 
to measurable time delays as a function of energy (for a review, see lAschwandenll2004l) : (1) 
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time-of-flight dispersion of free-streaming electrons, (2) magnetic trapping with the colli- 
sional precipitation of electrons, and (3) cooling of the thermal plasma. 



8.1 Time-of-Flight Delays 

The first type, the time-of-flight (TOF) delay, has a scaling of At(e) « e -1 / 2 and is caused 
by velocity differences of electrons that propagate from the coronal acceleration site to the 
chromospheric energy-loss region. The time differences ar e of order At w 10 — 100 ms fo r 
non-thermal electrons at energies E « 20- 100 keV (e.g.. lAschwanden et al.lfl995l.il 9961) . 
The measurement of such tiny time delays requires high photon statistics and high time 
resolution. Such data were provided by CG7?0/BATSE, which had 8 detectors, each with 
an effective collecting area of ~2000 cm 2 and oriented at different angles to the Sun so that 
detector saturation at high count rates was not a problem. (For comparison, the total effective 
collecting area of RHESSI's detectors is less than 100 cm 2 .) 

These studies of TOF delays have provided important evidence that electrons are ac- 
celerated in the corona, above the top of the hot flare loops observed in soft X-rays. The 
fine structure in the light curves of most, but not all, of the studied flare bursts showed 
energy-dependent time delays consistent with the free streaming of electrons to the foot- 
points of the flare loops from an origin somewhat more distan t than the half-length of t he 
loops ( Aschwanden et al. 1995; Aschwanden & Schwartz 1995; Aschwanden et al. 1996). 



8.2 Trapping Delays 

The second type, the trapping delay, is caused by magnetic mirroring of coronal electrons 
which precipitate toward the chromosphere after a collisional time scale At (e) °c £ 3 / 2 . This 
is observable f or time differences of At & I — 10s for non - thermal electrons at E w 20 — 
100 keV (e.g Jvilmer et al.ll 19821 ; lAschwanden et al.ll 19971) . For trapping delays the higher 
energy X-rays lag the lower energy X-rays, as opposed to time-of-flight delays where the 
hig her energy X-rays precede the lower energy X-rays. 

lAschwanden et all l l 19971) filtered variations on time scales ~1 s or less out of CGRO 
BATSE flare HXR light curves. They found time delays in the remaining gradually vary- 
ing component to be consistent with magnetic trapping and collisional precipitation of the 
particles. Trap plasma densities ~ 10 11 cm~ 3 were deduced. No evidence was found for a 
discontinuity in the delay time as a function of energy and, therefore, for second-step (two- 
stage) acceleration of electrons at energies < 200 keV. 



8.3 Thermal Delays 

The third type, the thermal delay, can be caused by the temperature depe ndence of cool- 
ing p r ocesses, such as by thermal conduction, T C (T) °= r~ 5 / 2 (e.g., Antiochos & Sturrock 



1978; 



1990; 



Culhane et al] 1 19941) . or by radiative cooling, X r (T) °= T 5 ^ (e.g., Fisher & Hawlev 
Cargill et alJ 1995l\ The observed physical parameters suggest that thermal conduc- 



tion dominates in flare loops at high temperatures as observed in soft X-ray wavelengths, 
while radiati ve cooling dominates in the later phase in postflare loops as observed in EUV 
wavelengths dAntiochos & Sturrocklll978l ; ICulhane et al.|[l994l ; lAschwanden & Alexander] 
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l200lh . When the temperature drops in the decay phase of flares, the heating rate can jus- 
tifiably be neglected and the conductive or radiative cooling rate dominate the temperature 
evolution. Before RHES SI, the cooling curve T(t ) in flare plasmas had been studied in only 
a few flares (e.g.. iMcTiernan et alj|l993l ICulhane et al.|[l994l ; lAschwanden & A lexander 
120011) . 

The high spectral resolution of RHESSI data is particularly suitable for any type of ther- 
mal modeling, because we can probe the thermal plasma from ~3 keV up to ~30 keV 
with a FWHM resoluti on of ~1 keV thanks to the cooled germanium detectors dlin et al] 
l200llSmith et al.l2002l) . This allows us to measure flare temperatures with more confidence. 
A statistical study of flare temperatures measured in the range of T w 7 — 20 MK indeed 
demonstrates some agreement between the values obtai ned from spectr al fitting of RHESSI 
data with those obtained from GOES flux ratios dBattaglia et al.ll2005t) . although RHESSI 
has a b ias for the high-temperature tail of the differential emissi on measure (DEM) distri- 
bution dAschwanden et al1l2008l : Ivaananen & Pohiolainenll2007h . Of course, we expect an 
agreement between the deduced emission-measure-weighted temperatures only when both 
instruments are sensitive to a temperature range that covers the flare DEM peak. 

A close relationship between the non-thermal and thermal time profiles was found early 
on, in the sense that the thermal emission often closely resembles the inte gral of the non- 
thermal emission, a relationship tha t is now known as the Neupert effect (Neu pertlll968l : 
lHudsonl[l99ll :l Dennis & Zarro 1993). This relationship is, however, strictly only expected 
for the asymptotic limit of very long cooling times, while a physically more accurate model 
would quantify this effect by a convolution of the non-thermal heating with a finite cooling 
time. The deconvolution of the e-folding cooling time in such a model has never been at- 
tempted statistically and as a function of energy or temperature. Theoretical discussions of 
the Neupert effect, inclu ding multiple energy release events, chromosp heric evaporation, and 
coolin g, can be found in lWarren & Antiochosl d2004l) . lLiu et alj d2010h . and lReeves & Moats! 

feoich . 

The cooling time at a given energy can be estimated from the decay time of a flare 
time profile. For instance, the decay t imes measured with G OES in soft X-rays were found 
to have a median of X^ ecay w 6 min dVeronig et ai]|2002d lbh. The observed cooling times 
have typically been found to be much longer t han predicted from classical conduction , 
but shorter than the radiative cooling time (e.g., McTierna n et al ] | 19931 : Ijiang et al.l[200l : 
iRavmond et al.ll2007h . This discrepancy could result from either continuous heating or sup- 
pression of conduction during the decay phase, or a combination of both. 

The Neupert effect was tested by correlating the soft X-ray peak flux with the (time- 
integrated) hard X-ray flux. A high correlation and time coincidence between the soft X-ray 
peak and hard X-ray e nd time was generally found, but a significant fraction of events also 
had a different timing dVeronig et alJl2002ch . A delay of 12 s was found in the soft X-ray 
flux time deriva tive with respect to the hard X-ray flux in SOL2003-11-13T05:01 (M1.6) 
dLiu et al] |2006l . also see Section PT2l >. Time delays such as this could be related to the hy- 
drodynamic flow time during chromospheric evaporation. Tests of the "theoretical Neupert 
effect," i.e., comparisons of the beam power supply of hard X-ray-emitting electrons and the 
thermal energy of evaporated plasma observed in soft X-rays, f ound it to strongly d epend 
on the low-energy cutoff to the non-thermal electron distribution dVeronig et all20 05). This 
provides another approach to deducing the energy at which the low-energy cutoff in the 
electron distribut ion occurs in indiv idual flares. The Neupert effect has also been studied in 
several flares bv lNindd200i,l2009h . who finds a high correlation between the hard X-ray 
flux and the time derivative of the thermal energy deduced from X-ray spectral fits (Ning 
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Fig. 8.1 Example of a multi-thermal spectrum with contributions from plasmas with temperatures of T = 
15, 20,..., 50 MK and a DEM distribution of dEM(T)/dT « T~ 4 . The individual thermal spectra and their 
sum are shown with thin linestyle, where the sum represents the observed spectrum. Note that the photons 
in the energy range e = 5.8 — 19.4 keV are dominated by temperatures of T=15-50 MK, which have a 
corresponding thermal energy that is about a factor of (4 + 1/2) = 4.5 lower than the corresponding photon 
energy (£,/, = 1.3 — 4.3 keV). The summed photon spectrum without the high-temperature cutoff approaches 
the power-law function 1(e) °= e~ 3 5 (dotted line) (from Aschwanden 2007). 

boost) and an anti-correlation between the hard X-ray spectr al index and the time rate of 
change of the UV flare area observed by TRACE dNing|l2009h . 



8.4 Multi-Thermal Delay Modeling with RHESSI 

Since major solar flares generally produce a large number of individual postflare loops, giv- 
ing the familiar appearance of loop arcades lined up along the flare ribbons, it is unavoidable 
that each loop is heated up and cools down at different times, so that a spatially integrated 
spectrum alw ays contains a multi-thermal differential emission measure distribution (cf. 
IWarrenll2006l) . The resulting multi-thermal bremsstrahlung spectrum (for photon energies 
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Fig. 8.2 X-ray light curves are shown for SOL2002-02-26T10:27 (C9.6), for energies of 10 keV to 30 keV 
in intervals of 1 keV, observed with RHESSI (left panels). The spectrum is decomposed into thermal and non- 
thermal components (top right panel) and the delay of the peaks at different energies is fitted with a thermal 
conduction cooling time model that has a scaling of T conl i(T) T~P (right bottom panel). The best fit shows 
a power index of j3 = 2.8, which is close to the theoretically expected value of j} = 5/2 (Equation l8.2t . The 
full delay of the thermal component is indicated with a thin curve (bottom left p anel), while the weighted 
(thermal+non-thermal) fit is indicated with a thick curve (from Aschwanden 20Q3). 



e) observed in soft X-rays (neglecting the Gaunt factor of order unity), 

W-t/ BffiQg^W , (8,, 

is then a function of a multi-thermal differential emission measure distribution dEM(T) = 
n 2 (T)dV. An example of a multi-thermal spectrum from a differential emission measure 
proportional to T~ 4 up to a maximum temperature of 50 MK is shown in Figure [87X1 

As discussed above, the initial cooling of the hot flare plasma (say at T> 10 MK) is 
generally dominated by conductive cooling (rather than by radiative cooling, which can 
dominate later after the plasma cools to EUV-emitting temperatures of T <2 MK). The 
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thermal conduction time has the following temperature dependence: 

r (t\ e ' h 3«Ar 2\ n e L 2 k B rr _ 5ll /ry 5/2 . 

see I Aschwandenl J2007t) for parameter definitions. Since the thermal bremsstrahlung at de- 
creasing photon energies is dominated by radiation from lower temperature flare plasma, the 
conductive cooling time is expected to become longer at lower temperatures (T C ond x T~ 5 I 2 ). 
Thus, the soft X-ray peak is always delayed with respect to the harder X-ray peaks, reflecting 
the c onductive cooling o f the flare loops. 

lAschwandenl d2007f) has measured and modeled this conductive cooling delay T coto /(e) 
for a comprehensive set of short-duration (< 10 min) flares observed by RHESSI. One ex- 
ample is shown in Figure [8^21 He finds that the cooling delay At expressed as a function of 
the photon energy e and photon spectral index y can be approximated by 



7 

Ar(e,r)«Tg- 



log 1 + — 



T g v(y-i)eo 



3/4 

(8.3) 



(where T g is the Gaussian width of the time profile peak) and yields a new diagnostic of 
the pr ocess of conductive cooling in multi-thermal flare plasmas. In a statistical study of 65 
flares ( Aschwanden 2007), 44 (68%) were well fit by the multi-thermal model with a best fit 
value for the exponent of jS = 2.7 ± 1.2, which is consistent with the theoretically expected 
value of jS =2.5 according to Equation 18.21 The conductive cooling time at 7b = 1 1.6 MK 
(£o = 1 keV) was found to range from 2 to 750 s, with a mean value of T c o = 40 s. 

We note that these timing data, as well as thick-target fits to the non-thermal part of spec- 
tra that reveal the evolution of the energy content in accelerat ed electrons, pr ovide additional 
constraints on models such as the multithread flare model of lWarrenl fiam . 



9 Hard X-ray spectral evolution in flares 

9.1 Observations of spectral evolution 

The non-thermal hard X-ray emission from solar flares, best observed in the 20 to 100 keV 
range, is highly variable. Often several emission spikes with durations ranging from seconds 
to minutes are observed. In larger events, sometimes a more slowly variable, long duration 
emission can be observed in the later phase of the flare. Hence, most flares start out with an 
impulsive phase, while some events, mostly large ones, show the presence of a late gradual 
phase in the hard X-ray time profile. 

While these two different behaviors can already be spotted by looking at lightcurves, 
they also are distinct in their spectral evolution. The impulsive spikes tend to be harder at 
the peak time, and softer both in the rise and decay phase. The spectrum starts soft, gets 
harder as the flux rises and softens again after the maximum of the emission. This pattern of 
the spectral evolution is thus called soft-hard-soft (SHS). On the other hand, in the gradual 
phase, the flux often slowly decreases, while the spectrum stays hard or gets even harder. 
This different kind of spectral e volution is called soft-hard-harder (SHH). 

Historicall y, both the SHS JParks & Wincklerlll969l : iKane & Andersonlfl970h and the 
SHH behavior dFrost & Dennislll97lh were observed in the early e ra of hard X-ray observa- 
tions of the Sun. Subsequent investigation confirmed both the SHS teendl977l ; lBrown & Loraril 
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Fig. 9.1 Time evolution of the spectral index y {upper curve, linear scale on right) and the flux normalization 
/35 (lower curve, logarithmic scale on left) of the non-thermal component in SOL2 002-1 1-09T13:23 (M4.9). 
Different emission spikes are shown in different colors (after Grids & Benz 2004). 



ll985l: lLin & Schwardll98llGanlll998l:lFletcher & Hudsodl200irHudson & Famikl| 2002 h 
and the SHH deliver et al.|l98d ; lKiplingerlll995l ; ISaldanha et al.l2008l : lGrigis & Benzl2 008) 
behavior. T he SHH behavior has been found to be correlated with proto n events in interplan- 
etary space dKiplingejll995l : ISaldanha et al.ll2008l : lGravson et al1l2009h . 

Evidence for hard-soft-hard (HSH) spectral evolution at energ ies above ~50 keV ha s 
been reported for multiple spikes in SOL2004-1 1-03T03:35 (Ml. 6) dShao & Huangll2009bl) . 
SHS behavior was observed at lower energies. This H SH behavior might be explained by 
albedo, which typically peaks around 30^-0 keV (see lKontar et alJfeoi ll) . but the authors 
corrected for albedo from isotropically emitted photons. A likely explanation is that the 
spikes overlie a harder, gradually varying component, possibly emission from trapped elec- 
trons (Section [82). 

While all these observations established the qualitative properties of the spectral evo- 
lution, a statistical analysis of the quantitative relation between the flux and spectral index 
had not been performed in the pre-RHESSI era. Here, we summarize RHESSI results inves- 
tigating quantitatively the spectral evolution of the nonthermal component o f the hard X-ray 
emission, as well a s the theoretical implications. More details can be found in lGrigis & Benj 
fe004l2005ll20o3) . 

To quantify the spectral evolution, a simple parameterization for the shape of the non- 
thermal spectrum is needed. Luckily, in solar flares the spectrum is well described by a 
power law in energy, which often steepens above 50 keV. Such a softening of the spectrum 
can be modeled by a broken power-law model. However, it is difficult to observe such a 
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downward bending at times of weak flux, because the high-energ y region of the spectrum 
is lost in the background. As a compromise. iGrigis & BenzU 2004) fitted the data to a single 
power-law function at all times. Although the single power law does not always provide 
a good fit to the spectra, it provides a characteristic spectral slope and ensures a uniform 
treatment of the spectra at different times. 

The two free parameters of the power-law model are the spectral index 7 and the power- 
law normalization 7 £o at the reference energy Eq. The reference energy £0 is arbitrary, but 
fixed, usually near the logarithmic mean of the covered energy range. In the RHESSI spectral 
analysis software, OSPEX, £0 = 50 keV by default. The time dependent spectrum is given 
by 

£ \ "TO 



7(8,0=^(01-) • (9-D 

A repre sentative sample of 24 solar flares of GOES magnitudes between Ml and XI was 
selected bv lGrigis & Benzl 1120041) . The spectral model (Equation [9TTJ, with the addition of 
an isothermal emission component at low energies, was fitted with a cadence of one RHESSI 
spin period (about 4 s). This delivered a sequence of measurements of the quantities 7^(0 
and y(0 for each of the 24 events, covering a total time of about 62 minutes of non-thermal 
hard X-ray emission. For these events, £0 = 35 keV was chosen, a meaningful energy which 
lies about in the middle of the range where the non-thermal emission is best observed in 
these M-class flares. 

An example of the measured time evolution of the spectral index 7 and the flux nor- 
malization 735 for the longer-lasting event of the set is shown in Figure 19.11 A correlation 
in time between the two curves can be readily seen. Single emission spikes are plotted in 
different colors, so that the soft-hard-soft evolution can be observed during each spike (with 
the exception of the late, more gradual phase, where the emission stays hard as the flux 
decays). 

As there is an anti-correlation in time between log 735 (f ) and y(t), a plot of one parameter 
as a function of the other, eliminating the time dependence, shows the relationship between 
them. Figure 19.21 shows plots of 7 vs. 735 for 3 events where there are only one or two 
emission peaks. The points in the longer uninterrupted rise or decay phase during each event 
are marked by plus symbols. A linear relationship between log 735 and 7 can be seen during 
each phase, although it can be different during rise and decay. 

On the other hand, a plot of all the 911 fitted model parameters for all the events show a 
large scatter, as shown in Figure [9731 The large scatter can be understood as originating from 
the superposition of data from a large numbers of different emission spikes, each featuring 
linear trends with different parameters. This plot does demonstrate, however, the tendency 
for flatter spectra to be associated with more intense flares. 

RHESSI observations of the gr adual phase of large solar flares 1 Grigis & Benz 200of) and 



its relation with proton events dSaldanha et al . 200 S lGravsonetal.ll2009h have shown that 
the hardening behavior is complex and cannot be characterized by a continuously increas- 
ing hardness during the event. Therefore the soft-hard-harder (SHH) denomination does not 
accurately reflect the observed spectral evolution. Rather, phases of hardening (o r even ap- 
proxi matively constant hardness) are often seen in larger events as the flux decays dKipl inger 
I1995I) . The start of the hardening phase can happen near the main peak of the flare, or later. 
The end of hardening can even be followed by new impulsive SHS peaks. The most re- 
cent s tatistical study of the correlation of SHH behavior with proton events dGravson et alj 
l2009h found that in a sample of 37 flares that were magnetically well-connected to Earth, 1 8 
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Fig. 9.2 Spectral index y vs. flux normalization 735 for three events, showing the linear dependence of single 
rise and decay phases of emission spikes on a log -linear scale. Dots ma rk results from individual spikes, while 
pluses mark the longer rise or decay phase (from Grieis & Benz 2004). 

showed SHH behavior and 12 of these produced solar energetic particle (SEP) events. None 
of the remaining 19 flares that did not show SHH behavior produced SEP events. 



9.2 Interpretation of spectral evolution 

Can we explain the soft-hard-soft spectral behavior theoretically? The problem here is that 
many effects contribute to the properties of the high-energy electron distribution whose 
bremsstrahlung hard X-rays are observed by RHESSI and similar instruments. We can iden- 
tify three main, closely related classes of physical processes that affect the distribution of 
the electrons and the spectrum of the X-ray photons they generate: (1) the acceleration of 
part of the thermal ambient plasma, (2) the escape from the acceleration region, and (3) the 
transport to the emitting region. The photon spectrum also depends on the properties of the 
bre msstrahlung emiss ion mechanism. 

lMilleretai]Jl996h proposed a stochastic acceleration mechanism where electrons are 
energized by small-amplitude turbulent fast-mode waves, called the transit-time damping 
model. They showed that their model could successfully account for the observed number 
and energy of electrons accelerated above 20 keV in subsecond spikes or energy release 
fragments in impulsive solar flares. However, they made no attempt to explain the observed 
hard X-ray spectra (which are softer than predicted by the transit-time damping model) and 
did not consider spectral evol ution. Furthermore, this approach does not account for particle 
escape. iGrigis & BerizU2006l) extended the model with the addition o f a term describing the 
escape of the particles from the acceleration region, as in the model of lPetrosian & Don aghv 
dl999h . To ensure conservation of particles, they also add a source term of cold particles 
coming into the accelerator (such as can be provided by a return current). 
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Fig. 9.3 Plot of the spectral index y versus the fitted non-therm al flux at 35 keV (give n in photons s 1 cm 2 
keV~'). All 91 1 data points from the 24 events are shown (from Grigis & Benz 2004). 



The stochastic nature of this acceleration model implies that the electrons undergo a 
diffusion process in energy space. Mathematically, the acceleration is described by the fol- 
lowing convective-diffusive equation: 

§7 = \ 7& i {DcOLL + Dj) f ] ~ Je i (Ac0LL + Aj) f ] ~ S{E } ' f + Q[E) ' (9 ' 2) 

where f(E) is the electron density distribution function, Df and Aj are, respectively, the 
diffusion and convection coefficients due to the interactions of the electrons with the accel- 
erating turbulent waves, Dcoll and Acoll are, respectively, the diffusion and convection co- 
efficients due to collisions with the ambient plasma, S(E) is the sink (escape) term, and Q(E) 
is the source (return current) term. The escape term is proportional to v(E) /t, where v(E) is 
the electron speed, and X is the escape time. The escape time can be energy-dependent, but 
for simplicity it is initially kept constant. The longer the escape time, the better the particles 
are trapped in the accelerator. The source term is in the form of a Maxwellian distribution 
of electrons with the same temperature as the ambient plasma. 

The coefficients Dj and At are proportional to the dimensionless acceleration parameter 

/acc = ^-^, (9.3) 

where Uj and Ub are, respectively, the energy densities of the turbulent waves and of the 
ambient magnetic field, (k) is the average wave vector, and Qn is the proton gyrofrequency. 
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Fig. 9.4 Accelerated electron density distributions with different values of the power-law index resulting 
from changes in 7 T = Ixcc ■ The dashed curve represents the ambient Maxwellian distribution. The two 
dotted lines indicate the energy range used for the co mputation of the pow er-law index 5 shown above each 
spectrum. Harder spectra have a larger 7 T value (from Grigis & Benz 200f|). 



Equation l9.2l can be solved numerically until an equilibrium state (df/dt = 0) is reached. 
The equilibrium electron spectra from the model are controlled by two parameters: the ac- 
celeration parameter 7acc described above and the escape time T. Above 10-20 keV, the 
collision and source terms in Equation ( 19.2b can be neglected, since they apply to the am- 
bient Maxwellian, and thus the equilibrium spectra depend to a first approximation only on 
the product 7 T = /acc ■ 

Figure |9~4l shows the equilibrium electron spectra for different values of l t = I acc • T. 
As I t increases, the spectrum gets harder and harder. To explain the soft-hard-soft effect, 
either the acceleration or the trapping efficiency (or both) must increase until the peak time, 
and then decrease again. We note that this model does not include magnetic trapping (other 
than in the magnetic turb ulence itself), which can alt er the computed electron spectra and 
their time evolution (e.g., iMetcalf & Alexandejfl999T) . 

To see whether this pr oduces the linear relat ion between the spectral index and the log 
of the flux normalization, iGrigis & Beiizl J2006h computed the hard X-ray emission from 
these model electron spectra. Since these are equilibrium spectra, thin-target emission was 
computed. They then plotted the spectral index vs. the flux normalization of the resulting 
photon spectra. Since the spectra are not power-law, but bend down, they fit a power-law 
model to the model photon spectrum in a similar range as the one used for the observations. 

Figure 19.51 shows the computed values for the spectral indices and flux normalizations 
for both the electron and the photon spectrum from the model. The results show that there 
is indeed a linear relation between the spectral index and the log of the flux normalization. 

An alternative mechanism that could be responsible for soft-hard-soft spectral evolution 
is return current losses as the electrons propagat e to and within the thick-target footpoints 
of the flare loop l lZharkova & Gordovskvvl [2006). The highest electron energy to which re- 
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Fig. 9.5 Model results for the spectral index and flux normalization for electrons and photons. The dashed 
line is the best straight-line fit to the model results (in the range of spectral i ndices from 2 to 8 for the electrons, 
and 3 to 9 for the photons), corresponding to a pivot-point behavior (from Grigis & Benz 2006). 



turn current losses are significant is proportional to the return current electric field strength, 
which is in turn proportional to the electron beam flux density (see Section|5]l. Therefore, as 
the electron flux density increases and then decreases, the low-energy part of the X-ray spec- 
trum flattens to higher and then lower energies as the return current electric field strength 
increases and then decreases. The net effect is SHS spectral evolution below the maximum 
energy for which return current losses are significant during the flare. The observation of 
SHS behavior in coronal X-ray sources, however, indicates that this spectral evolution is a 
property of the accel eration mechanism rathe r than a consequence of energy losses during 
electron propagation dBattaglia & BenzfeOOrj . see Section [l0.3b . 

Are there two stages of electron acceleration, one responsible for the impulsive phase 
and one for the gr adual phase? RHESS I spectroscopy and imaging of a set of 5 flares with 
hardening phases ( iGrigis & Bend l2008) showed that there is no discontinuity in the motion 
of footpoints at the onset of hardening and no clear separation between the impulsive and 
the gradual phase: the former seems to smoothly merge into the latter. This supports the 
view that the same acceleration mechanism changes gradually in the later phase of the flare, 
rather than a two stage acceleration theory. The hardening phase may in fact be caused by 
an increase in the efficiency of trapping of the electrons above 100 keV. 

The underlying cause of the S HS spectral evolution has b een a ddressed in terms of th e 
stochastic acceleration model by Bvk ov & Fleishman! d2009h and lliu & Fletcher (2009). 
Bykov and Fleishman consider acceleration in strong, long-wavelength MHD turbulence, 
taking into account the effect of the accelerated particles on the turbulence. They argue 
that the electron spectrum flattens during the linear acceleration phase, while the spectrum 
steepens during the nonlinear phase when damping of the turbulence because of the parti- 
cle acceleration is important, giving SHS spectral evolution. They argue that SHH evolu- 
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tion will be observed when the injection of particles into the acceleration region is strong. 
Liu & Fletcher also argue that the SHS evolution results from dependence of the electron 
distribution power-law index on the level of turbulence as it increases and subsequently 
decreases. They attribute changes in the SHS correlation during a flare to changes in the 
background plasma, likely due to chromospheric evaporation. 

We note that simple direct-current (DC) electric field acceleration of electrons out of the 
thermal plasma can produce the SHS spectral evolution. The flux of accelerated electrons 
and the maximum energy to which electrons are accelerated and, therefore, the high-energy 
cutoff to the electron distribution, increa se and decrease together as the electric field strength 
increases and decreases ( lHolmarj[l985h . The X-ray spectru m is steeper a t energies within 
one to two orders of magnitude below the high-energy cutoff (Holman 2003). In large flares, 
however, where the X-ray spectrum is observed to continue to MeV energies or higher, there 
is no evidence for a high-energy cutoff in the appropriate energy range. Therefore, at least for 
large flares with spectra extending to high energies, a simple DC electric field acceleration 
model does not appear to be appropriate. 



10 The connection between footpoint and coronal hard X-ray sources 

Hard X-ray (HXR) sourc es at both footpoints of a coronal loop structure have been observed 



since Hovng et al. Ul98ll) . As reviewed in Sections [T] & [2j they are understood to be thick- 
target bremsstrahlung emission produced by precipitating electrons, accelerat ed somewhere 
in or a bove the loop. A third HXR so urce situated above t he looptop (see iKrucker et al] 
2008a, for a review) was first noted by iMasuda et al.l ll 19941) in Yohkoh observations. The 
nature of this coronal HXR source has remained uncertain, but in simple solar flare models 
with reconnection and particle acceleration in the corona, we expect some relation between 
coronal HXR sources and footpoints. RHESSI has enabled us to study events featuring coro- 
nal HXR sources and footpoints simultaneously. By studying the behavior of the sources in 
time and the relations between them, we can address questions such as: Are both coronal 
and footpoint emissions caused by the same electron population? How is such an electron 
beam modified in the loop (collisions, return currents, trapping, etc.)? Is SHS behavior (Sec- 
tion ^. U a transport effect produced by collisions or return currents, or is it a feature imposed 
by the acceleration mechanism? 



10.1 RHESSI imaging spectroscopy 

RHESSI has provided the possibility of obtaining simultaneous, high-resolution imaged 
spectra at different locations on the Sun. One can therefore study each source separately 
in events with several contemporaneous HXR sources. The high spectral resolution has al- 
lowed a reliable differentiation between thermal and non-thermal emission to be made in 
many flares. Furthermore, RHESSFs imaging spectroscopy has allowed differences in indi- 
vidual flare source spectra and their evolution to be studied in considerable detail. 

Imaged spectra and the relative timin g of sources in three fla r es, including the limb flare 
SOL2002-02-20T11:07, were studied by krucker & Linl ( l2002h . ISui eTal] ( l2002h analyzed 
and modeled the tw o footpoint sources and a high, above-the-looptop hard X-ray source ob- 
served in this flare. lEmslie et al.l d2003l) analyzed SOL2002-07-23T00:35 (X4.8) flare with 
four HXR sources observed by RHESSI. They found a coronal source with a strong thermal 
component, but the non-thermal component could not be studied due to severe pulse pile-up. 
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iBattaelia & Beiizl ((2006) studied five M-class events. Due to the smaller pile-up amount in 
those events, studying the non-thermal coronal emission was possible. The results of these 
studies are summarized below. 



RHESSI 34-38 keV 13-Jul-2005 14:14:00.000 UT 




10 100 10 100 

Energy (keV) Energy (keV) 

Fig. 10.1 Top left Composite CLEAN image of a RHESSI event with three hard X-ray sources. The footpoints 
(labeled 1 & 2) are visible on the solar disc in an image made at 34-38 keV. The position of the coronal source 
(labeled 3) high above the limb is indicated by the 50 and 80% white contours taken from a 10-12 keV image. 
Plots 1-3 show spectra and normalized residuals over the fitted energy range for the north footpoint (1), south 
footpoint (2), and coronal source (3) (after lBattaglia & Benzl200U) . 



10.2 Relation between coronal and footpoint sources 

The quantitative relations between the footpoints and the coronal source and between the 
two footpoints can give information about the physical mechanisms at work in a solar flare. 
Simple models envision a beam of accelerated electrons encountering a low-density region 
in the corona, leading to thin-target bremsstrahlung. When the same electron beam reaches 
the chromosphere, the particles are fully stopped in the dense material, producing thick- 
target emission. Assuming an electron power-law distribution for the electron energy E of 
the form 

^(E)=AE~ S (10.1) 
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producing thin-target bremsstrahlung in the coronal source, the observed photon spectrum 
has spectral index jthin = 8+1 (Equation [2T0}. Reaching the chromosphere, the accelerated 
electrons will be fully stopped, producing thick-target bremsstrahlung with a photon spectral 
index ythick = 8 — 1 (Equation l2.12b . In such a simple scenario one would therefore expect 
a difference in the photon spectral index Ythin ~ Ythick = 2 between the coronal source and 
the footpoints. Further, the two footpoints should be of equal hardness and intensity if one 
assumes a symmetric loop and symmetric injection of particles into the legs of the loop. 

10.2.1 Observed difference between coronal and footpoint spectral indices 

A sample of flares o bserved with Yohkoh to have coronal HXR sources was studied by 
Petrosian et ail J2002h . They found that the spectral index of the coronal sources was, on 
the average, steeper by 1 than the spectral indices of the footpoint sources. ISui et alJ ( Eoo3) 
also found a spectral index difference of 1 for SOL2002-02-20T 11:07 (C7.5) observed with 
RHESSI. 

iBattadia & Bend ll200fih found that the coronal source was softer than both footpoints 
for all of their five events in nearly all analyzed time bins. Figure [ToTI (top left) shows an 
image of SOL2005-07-13T 14:49 (M5.0) in the 34-38 keV energy band. The two footpoints 
are visible, as well as the 50 and 80% contours of the coronal source taken from a 10- 

12 keV image. Spectra and spectral fits are shown for the two footpoints and the coronal 
source. The steepness of the coronal source spectrum (number 3 in the figure) relative to 
the spectra from the footpoints is apparent. However, the quantitative difference between 
the values of the spectral index obtained for the coronal source and the footpoints often 
differed significantly from 2. For the five flares analyzed, the smallest mean difference in 
the spectral indices, averaged over time, was 0.59±0.24. The maximum mean difference, 
averaged over time, was 3.68±0.14. These clearly contradict the theoretical expectation 
summarized above. Simple thin-thick target scenarios do not seem to work in most cases 
and additional effects need to be considered. 

Evi dence for two populations of coronal source non-thermal spectra was found bv lShao & Huand 
(2009a). They compare coronal and footpoint spectral indices at 28 hard X-ray peaks from 

13 single-loop flares observed by RHESSI. The spectral index in the coronal sources was 
determined from an isothermal plus power-law fit below 30 keV, while the footpoint spec- 
tral indices were determined from a power-law fit at 30-60 keV photon energies. They argue 
that the coronal spectra can be divided into two groups. One, for which the coronal spectral 
index is greater than 5, is well correlated with the footpoint spectral index, and the difference 
in the indices ranges from 2-4. For the other, where the spectral indices are anticorrelated, 
the coronal spectral index is less than 5, and the difference in the indices ranges from 0-2. 
For the group of anticorrelated spectral indices, the coronal spectral index is correlated with 
the photon flux, while the footpoint spectral index is anticorrelated with the photon flux for 
both groups. These are intriguing results if confirmed by future studies. 

10.2.2 Differences between footpoints 

No significant differenc e was found in the spec tral i ndices for t he two footpoints in SOL2 002- 
02-20Tll:49 (C7.5) by iKrucker d2002h and lSui et all fc002b . IPiana et al l d2007h in- 

verted count visibility spectra for this flare to obtain mean electron flux distributions for 
the footpoints. They found the mean electron flux distribution function at the northern foot- 
point to be somewhat steeper (AS « 0.8) than that derived for the southern footpoint. They 
also found the distribution function for the region between the footpoints (not the coronal 
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source studied by Sui et al.) to be steeper than the footpoint distribution functions (A8 « 1.6 
rel ative to the southern fo otpoint) and to substantially steepen at energies above ~60 keV. 

lKrucker&Linl(l2002h found that, when a connection between footpoints could be deter- 
mined, the footpoints brightened simultaneously (to within the ~1 s time resolution of the 
observations) and had similar spectra. 

Differences of 0.3 - 0.4 betwe en the spectra l indic es of two footpoints in SOL2002-07- 
23T00:35 (X4.8) were reporte d bylEmslie et al.l d2003l) . 

For the flares analyzed bv lBattaglia & B enz (2006), a significant difference was found 
in only one out of five events. For all other events, the mean difference in 7f p was zero within 
the statistical uncertainty. 

Different spectra at the two footpoints imply an asymmetric loop. Such an asymmetry 
can result, for example, from different column densities or different beam fluxes and cor- 
responding return current energy losses in the leg s of the loop. It could also result from 
asymmetric magnetic trapping within the lo op (e.g., I Alexander & Metcali 2002). In a study 
of 53 flares showing two HXR footpoints, ISaint-Hilaire et al] d200St) found that footpoint 
asymmetry was greatest at the time of peak HXR flux and the difference in the footpoint 
spectral indices A y rarely exceeded 0.6. In most cases they found the footpoint asymmetry 
to be inconsistent with different colum n densities in the two legs of the loops. 

In SOL2003-10-29T20:49 (XlO.O) lLiu etaD J2009ah found that the brighter HXR foot- 
point was marginally, but consistently harder than the dimmer footpoint by A y = 0.15 ± 
0.13. They concluded that neither asymmetric magnetic mirroring nor asymmetric column 
density alone can explain the full time evolution of the footpoint HXR fluxes and spectral 
indices. However, a self-consistent explanation might be obtained by considering these two 
effects together and/or in combination with one or more additional transport effects, such as 
nonuniform target ionization, relativistic beaming, and return current losses. 



10.3 Spectral evolution in coronal sources 

Previous observations of SHS spectral evolution (see Section [97Tb were made with full-Sun 
spect ra which, except for ove r-the-limb events, are typically dominated by footpoint emis- 
sion. iBattaelia & Bend d2006l) . in their imaging spectroscopy study, found that the coronal 
source itself shows SHS evolution. This is illustrated in Figure 110.21 This finding implies 
that SHS is not caused by transport effects within the flare l oop, b ut is rather a property 
of the acceleration mechanism itself. Indeed, iGrigis & Berizl (2006) showed that SHS can 
be reproduced for electron spectra in a transit-time-damping, stochastic-acceleration model 
(Section[972l. 



10.4 Interpretation of the connection between footpoints and the coronal source 

In the above account, emphasis was given to the difference in the spectral index between the 
coronal source and footpoints. Assuming a thin target in the corona and a thick target at the 
footpoints, one would expect a difference of two. However, whether the coronal source acts 
as thin- or thick-targe t depends on the energy o f the accelerated electrons and the column 
density in the corona. IVeronig & Brownl ( |2004|) . for example, found coronal sources with 
column densities high enough to act as th i ck tar gets for electrons with energies up to 60 keV. 

As early as 1976. iMelrose & Brownl d 19761) showed that magnetic trapping with colli- 
sional scattering of electrons out of the trap can lead to a thick-target coronal source. The 
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Fig. 10.2 Top: GOES 1-8 A light curve of SOL2003-10-24T02:54 (M7.6). Middle: RHESSI 25-50 and 50- 
100 keV light curves near the peak of the GOES flare. Bottom: time evolution of fitted coronal source flux 
at 35 keV (F35, * symbols, left log scal e) and spectral index (/, + symbols, right log scale) displaying SHS 
evolution fafter lBattaglia & Beni^wS) . 



coronal source transitions through a thin-thick period, with the time scale for this transition 
depending on the electron energy and the plasma den sity in the trap. The trapping e ssentially 
increases the effective column density in the corona. iMetcalf & Alexande rl dl999h analyzed 
six flares with coronal sources observed by Yohkoh and found that three of the six flares 
showed properties consistent with trapping. 

A simple 1-D model that described the cor onal emission as intermediate thin- thick, de- 
pending on electron energy, was developed bv I Wheatland & Melrose! dl995T) . In this model 
a high-density region (>10 12 cm~ 3 ) is hypothesized to be present at or above the top of the 
flare loop. The model makes predi ctions for the sh ape of the coronal and footpoint spec- 
tra and the relations between them. Fletcte] dl995h obtained Monte Carlo solutions to the 
Fokker-Planck equation to show that, with the inclusion of high electron pitch angles and 
collisional scattering, a compact coronal X-ray source is produced at the top of a loop with a 
constant coronal density ~3 x 10 cm - 3 . lHolmanl ( fl99r3) showed that, even in the simple 1- 
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D model, a compact coronal source is produced when electrons are injected into a loop with 
a constant coronal density ~2 x 10 11 cirT 3 (see hesperia.gsfc.nasa.gov/sftheory/loop.htm). 
A compact coronal HXR so urce can also be produced i f there is a compact magnetic trap at 
or above the top of the loop. iFletcher & Martens! ( fl998h showed that, with such a trap, a sig- 
nificant coronal X-ray source can be produced at plasma densities as low as ~4 x 10 9 cm~ 3 . 
IPetrosian & Donaehvl dl999h showed that the coronal HXR source can be a consequence 
of acceleration and trapping by turbulence or plasma waves. In their stochastic accelera- 
tion model, the difference between the coronal and footpoint spectra is explained by the 
energy-dependent time scale for electrons to escape the acceleration region. 
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Fig. 10.3 Left: spectra fo r coronal source (red) and footpoints (blue) according to the model of 
Wheatland & Melrose 1 1995). The spatially integrated spectrum is shown in violet. Right: observed RHESSI 
spectra for the event SOL2()03-10-24T02:54 (M7.6). Isothermal and power-law fits to the coronal (crosses) 
and footpoint (dots) spectra are shown. T he vertical line indicat es the predicted critical energy for the transi- 
tion between thin and thick target (after Battaglia & Benz 2007). 



The left panel of Figure [ITHl illustrates the model of lWheatland & Melrosel Jl995h . The 
spatially integrated spectrum (violet) is the power-law spectrum (thick-target, Ymck = 5—1) 
expected for a single-power-law electron distribution with no low- or high-energy cutoffs 
and no thermal component. For £ « £ c = \/2KN (see Equation 12.4b , the spectrum is dom- 
inated by thick-target radiation from the coronal source (red). There is a low-energy cutoff 
in the electron distribution at the footpoints at E ps y/2KN because of the energy losses in 
the coronal source. The spectrum is dominated by thick-target radiation from the footpoints 
(blue) where £ E c . It is in this regime that the radiation from the coronal source is thin- 
target and the spectral index of the coronal source is steeper by 2 than that of the footpoints. 
These spectra are characteristic of all the models reviewed above. 

ISui et alJ d2002h compared the RHESSI observations of SOL2002-02-20T11:07 (C7.5) 
to a model with a consta nt-coronal-density loop and no magnetic trapping. They used a finite 
difference method (e.g jMcTiernan & Petrosianll990l : lHolman et alj2002i) to obtain steady- 
state solutions to the Fokker-Planck equation with collisional scattering and energy losses. 
Model images were convolved with the RHESSI response to produce simulated RHESSI 
observations for direct comparison with the SOL2002-02-20T1 1:07 flare images and imaged 
spectra. They found that, after obtaining a power-law model spectrum with an index of y = 3 
that agreed with the observed footpoint spectra, the effective spectral index of the coronal 
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source from the model (y = 4.7) was significantly steeper than that obtained for the flare 

(y=i). 

iBattaglia & Benzl d2007h compared the model of Iwheatland & Melrose] lQj>95) to the 
results of their study of five flares observed by RHESSI. The right panel of Figure 110.31 
shows observed spectra and spectral fits for one particular event. The observed spectra were 
dominated by thermal coronal emission at low energies. Therefore, not all of the model 
predictions could be tested. However, the observed relations between the spectra did not 
agree with the predictions of the model. For the flare in Figure [To31 for example, the differ- 
ence between the coronal source and footpoint spectral indices at the higher photon energies 
is 3.8 ±0.1, not 2. Also, an estimate of the column density in the coronal source gives 
\/2KN w 10-15 keV, while the interse c tion o f the coronal and footpoint spectra is found 
to be at e w 23 keV. Battaglia & Benz (2008) have found that this large difference in the 
spectral indices is consistent with spectral hardening caused by return current losses (see 
Section[5]>. 



11 Identification of electron acceleration sites from radio observations 

While energetic electrons excite hard X-ray emission during their precipitation into the 
dense layers of the solar atmosphere, they can also excite decimeter and meter wave ra- 
dio emission during propagation and trapping in magnetic field structures in the dilute solar 
corona. The radio emission pattern in dynamic spectrograms can give information about the 
electron acceleration process, the locations of injection of electrons in the corona, and the 
properties of the coronal magnetoplasma structures. 




Fig. 11.1 SOL2003-10-28T11:10 (X17.2). Left, bottom: fsee.Wa rmuth et all2007l) 200-400 MHz radio spec- 
trum (Astrophysical Institute Potsdam) showing the signature of the outflow termination sh ock (TS, starting a t 
1 1 :02:47 UT). Left, top: INTEGRAL count rates at 150 keV and 7.5-10 MeV. Right: (after lAurass et alj2007l) : 
radio source positions (Nancay Radio Heliograph, 327 MHz) overlaid on a SOHO-EYI image (11:47 UT 
195 A). The bright areas are EUV flare ribbons in AR10486. RHESSI HXR centroids are shown as "+". The 
integration time intervals are: for the TS source SW of AR10486 1 1:02:45-1 1:03:15 UT, for the continuum 
source CONT N of AR10486 11:13-11:17 UT, respectively (see also Figure \\L2\ . The radio contours are at 
50, 70, and 99.5% of the peak flux value. 

Here we take as an example SOL2003-10-28T1 1:10 (X17.2) (shown in Figure [TTTT i. 
Different acceleration sites can be discriminated during the impulsive and the gradual flare 
phases. Radio spectral data from the Astrophysical Institute Potsdam (AIP; iMann et alj 
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Il992h . imaging data from the Nancay Radio Heliograph (NRH. lKerdraon & Delouislll997l) . 
and hard X-ray (RHESSI, INTEGRAL) data were combined in the analysis of this event. 
The conclusion was reached that a nondrifting, high-frequency type II radio burst signa- 
ture in the radio spectrum coincided with a powerful electron acceleration stage. Simul- 
taneously with the nondrifting type II signature, highly relativistic (>10 MeV) electrons 
were observed in the impulsive phase of the flare (Figure fTl.il upper left). The radio spec- 
trum sugge sts that this can be du e to acceleration at the reconnection outflow t ermina- 
tion shock I Aurass & Manr] 12004 ), as predicted for a classical two-ribbon flare ( iForbesI 
1 19861 . iTsuneta & Naitdll998l . lAurass et a"i1l2002l) . The radio source site is observed about 
210 Mm to the SW of the flaring active region (TS in Figure fTTTTl right). In this direc- 
tion, TRACE and SO/ZO/LASCClj C2 images reveal dynamically evolving magnetoplasma 
structures in an erupting arcade ( I Aurass et al .11 20061) . For realistic parameters derived from 
these obser vations (the geome try, density, temperature, and low magnetic field values of 
~5 Gauss). iMann et aD 1 2006) have found that a fully relativistic treatment of acceleration 
at the fast-mode out flow shock can explain the observed fluxes of energetic particles (see 
IZharkovaetai]|201lh . 





Fig. 11.2 Timing of the source CONT in Figure fTTTTl the NRH 327 MHz flux curve (asterisks) versus the 
GOES 0.5-4 A flux curve (solid line, partly off-scale). Inset: SOHO/EYT image showing the radio source 
centroid (white asterisk) and RHESSI HXR centroids as in Figure 111.11 Thick bar: the start time of GeV- 
energy proton injection in space (after Aurass et al. 20Q7|). 



In the main flare phase of the same event, an additional radio source (CONT in Fig- 
ure QTTTJl was found, lasting for ~10 min, indicating the presence of another acceleration 
site. No X-ray, EUV, or Ha emission was observed at the location of this radio source. Fig- 
ure 111.21 gives the timing and the source position with respect to the flaring active region. 
CONT is a m-dm-continuum source with fiber burst fine structure. Fiber bursts are excited 
by whistler waves propagating along field lines of the coronal magnetic field. As marked by 
a bold bar in the Fi gure, the time of the CONT emission is also the start time of GeV proton 
injection in space. lAurass et al. I d2006h have shown that this source site is not far from an 
open field (particle escape) region in the potential coronal magnetic field. The source briefly 
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flash es up already in the early impulsive ph a se. Based on a new method of fiber burst analy- 
sis dAurass et al.l2005l : lRausche et alj2007h . lAurass et alJ d2007h argue that this source most 
likely indicates acceleration at a contact between separatrix surfaces of different magnetic 

flux systems. 

Ra dio observations of flares and their implications are further addressed in I White et alJ 



12 Discussion and Conclusions 

12.1 Implications of X-ray observations for the collisional thick-target model 

As discussed in Section [2] the core assumption of the collisional thick-target model is that 
the spatially integrated hard X-ray emission from non-thermal electrons is bremsstrahlung 
(free-free radiation) from electrons that lose all their suprathermal energy through colli- 
sional losses in the ambient plasma as they simultaneously radiate the hard X-rays. "Simul- 
taneously" means within the observational integration time. This implies that all electrons 
that contribute significantly to the observed radiation reach a plasma dense enough or, more 
precisely, traverse a high enough column density for all of their suprathermal energy above 
the observed photon energies to be collisionally lost to the ambient plasma within the inte- 
gration time. For typical >1 s integration times, these conditions are met when the electrons 
stream downward from the corona into the increasingly dense plasma of the solar transition 
region and chromosphere. 

Since the thick-target model is often implicit in our interpretation of the hard X-ray 
emission from flares, it is important to keep the underlying assumptions in mind and test the 
model while at the same time applying it to flare observations. We have discussed above sev- 
eral physical processes that, if significant, change the conclusions of the simple collisional 
thick-target model regarding the electron distribution produced in the acceleration region. 
These processes occur in either the thick-target region itself, or during the propagation of the 
electrons from the acceleration region to the thick-target region. Only with the high spec- 
tral resolution and imaging of RHESSI has it become possible to observationally address 
these processes. Even with the RHESSI observations, however, it is difficult to conclusively 
determine the importance of each process. 

A physical pro cess that distorts the emitted X-ray spectrum is albedo (Section [3 .41 and 
lKontaretal.ll201ll) . Fortunately, the albedo contribution to the X-ray spectrum can be easily 
corrected on the assumption that the X-ray photons are isotropically emitted. This correction 
is available in the RHESSI spectral analysis software. If the photons are significantly beamed 
downward, however, the distortion of the spectrum can be substantially greater than that 
from isotropically emitted photons. An anisotropic photon distribution results from emit- 
ting electrons with an anisotropic pitch-angle distribution. The degree of anisotropy of the 
electron pitch-angle distribution also quantitatively affects conclusions from the thick-target 
model concerning the acceleration process. Therefore, it is important to better determine the 
pitch-angle distribut ion of the emitting electrons and the contribution of albedo to the hard 
X-ray spectrum (see lKontar et al 

The simple collisional thick-target model assumes that the target plasma is fully ionized. 
We have seen, however, that a nonuniformly ionized target region can produce an upward 
kink, or "chicane," in an otherwise power-law X-ray spectrum (Section [4}. This spectral 
shift can provide a valuable diagnostic of the ionization state of the target plasma and its 
evolution. It is likely, however, that the power-law spectrum below the chicane is hidden by 
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thermal radiation. The chicane is then observed only as a downward break in the spectrum at 
energies above those dominated by the thermal emission. The upper limit on the magnitude 
of the break provides a method for ruling out nonuniform ionization as the sole cause of large 
spectral breaks. To further distinguish this break from spectral breaks with other causes, it 
is important to better determine the degree of ionization as a function of column density at 
the thick-target footpoints. 

Return-current energy losses can also produce a downward break in the X-ray spectrum 
(Section |5J- The break energy depends on both the thermal structure of the plasma in the 
flare loop and on the non-thermal electron flux density distribution. These spectral modifi- 
cations and their evolution throughout flares provide an important test for the presence of 
initially un-neutralized electron beams and the return currents they must drive to neutralize 
them. Although RHESSI observations provide substantial information about the structure 
and evolution of flare spectra, only a lower limit on the electron flux density can usually 
be determined. Observations and analysis sufficiently accurate and comprehensive to verify 
the presence of return current energy losses as the cause of a spectral break are yet to be 
obtained. On the other hand, significant evidence exists (Section \53§ indicating that return 
current losses do have an impact on flare hard X-ray emission. 

A thorough comparison of measured flare spectra with theoretical spectra computed 
from models incorporating collisional and return current energy losses (including their ef- 
fect on the angular distribution of the non-thermal electrons), as well as nonuniform target 
ionization and albedo, is still needed. Spectral fitting alone, however, is not likely to distin- 
guish the importance of these different mechanisms. Comparison of the time evolution of 
the spectra, as well as of the spatial structure of the X-ray emission, with expectations would 
certainly enhance the possibility of success for such an endeavor. 

The analysis of the evolution of X-ray source positions and sizes with photon energy 
and time provides another important test of the collisional thick-target model (Section [TJ. 
For these flares that show non-thermal source evolution in the corona and upper transition 
region, the source position and size are sensitive to the energy losses experienced by the 
non-thermal electrons. They are, in fact, sensitive to the very assumption that the sources are 
produced by electrons as they stream downward from an acceleration region higher in the 
corona. Further studies of the evolution of these coronal X-ray sources should substantially 
clarify the applicability of the collisional thick-target model. 

For completeness, we note that under some circumstances other radiation mechanisms 
may significantly contribute to the X-ray emission from non-thermal electrons. The possibil- 
ity that recombination (f ree-bound) radiation from the non-ther mal electrons is som etimes 
important is discussed in brown & Mallikl J2008l . l2009h (also see lKontar et al.ll201lh . How- 
ever, the contribution of non-thermal free-bound radiation has recently been found to b e less 
significant than originally estimated ( Brown et al koiCh . lMacKinnon & Mallikl j201Ch have 
concluded that inverse Compton radiation may significantly contribute to the X-ray/y-ray 
emission from low-density coronal sources. 

Another testable aspect of the collisional thick-target model is the heating of the flare 
plasma by the non-thermal electrons. If the flare plasma is primarily heated by these elec- 
trons and the thick-target region is primarily in the chromosphere and lower transition re- 
gion, heating originating in the footpoints and expanding into the rest of the flare loop 
through "chromospheric evaporation" should be observed. On the other hand, if the loop is 
dense enough for the thick-target region to extend into the corona or if return-current heating 
is important, localized coronal heating and different ion abundances should be observed. 

It has generally been difficult to establish a clear connection between the location and 
evolution of X-ray sources produced by non-thermal electrons and by thermal plasma at 
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different temperatures. This is largely because of a lack of high-cadence EUV images cov- 
ering a broad range of coronal and transition region temperatures prior to the launch of 
SDO. Future studies of the coevolution of non-thermal X-ray sources and thermal sources in 
flares will be important in determining the extent to which heating mechanisms other than 
collisional heating by non-thermal electrons is significant. 

Predicting the expected evolution of the heated plasma is hampered by insufficient 
knowledge of the dominant heat transport mechanisms. We have seen evidence that many 
flares cool by classical thermal conduction or radiation once the heating has subsided (Sec- 
tion l8.4b . but this is not likely to be the dominant transport mechanism during rapid heating. 
Nevertheless, the spatial evolution of flare X-ray sources has so far been found to be con- 
sistent with chromospheric evaporation (Section 17.2b . Also, the Neupert effect, observed 
in most flares, and Doppler-shift measurements qualitatively support the thick-target model 
(Section 18.3b , but these do not rule out the possibility of other heating mechanisms tem- 
porally correlated with the electron beam collisional heating. As discussed in Section [3] 
substantial progress has been made in deducing the energy flux (total power) carried by 
non-thermal electrons, but we usually can deduce only a lower limit to this energy flux. Con- 
tinuing studies of flares similar to SOL2002-04-15T03:55 (Ml. 2) and the initially cooler, 
early-impulsive flares (Section [331 > may provide a better handle on this energy flux for com- 
parison with thermal e volution. The therma l properties, energetics, and evolution of flares 
are discussed further in lFletcheretaD ( l201lh . 



12.2 Implications of X-ray observations for electron acceleration mechanisms and flare 
models 

In Section \W[ we addressed the X-ray spectra of hard X-ray sources sometimes observed 
above the top of the hot loops or arcades of loops observed in flares. We reviewed results 
indicating that the spectra are qualitatively, but not quantitatively consistent with expecta- 
tions for electrons passing through a thin-target or quasi-thick-target region on their way to 
the thick-target footpoints of the flare loops. The apparent failure of these relatively simple 
models is probably a manifestation of the more complex above-the-looptop X-ray source 
structure revealed by RHESSI observations. 

Before RHESSI, time-of-flight delays in hard X-ray timing indicated that electrons were 
accelerated in a region somewhat above the looptops of the hot flare loops in most flares 
(Section 18.1b . Also, cusps were observed at the top of flare loops by Yohkoh (e.g., Sec- 
tionllOb. indicating a magnetic connection to the region above the hot loops. 

RHESSI images have revealed flares with double coronal sources, one at or just above 
the top of the hot loops and the other at a higher altitude above the lower source. The 
centroid of the lower source is higher in altitude at higher X-ray energies, while the cen- 
troid of the upper source is lower in altitude at high er X-ray energies, indicating that en- 
ergy release occurred between these coronal sources dSui & Holmanll2003l ; ISui et alj|2004l ; 
Liu et afl2008tl2009bh . In one flare, the upper source accelerated outward to the speed of a 
subs equent coronal ma ss ejection. The white-light coronagraph on the Solar Maximum Mis- 
sion ([Webb et alj|2003l) . the Large Angle and Spectrometric Coronagraph (LASCO) and th e 
Ultraviol et Coronagraph S pectrometer ( UVCS) on SOHO dKo et alKollLin et alj|2003) . 
RHESSI dSui et aljl2005bh . and TRACE dSui et alJ l2006a) have all provided direct evidence 
for the presence of an extended, vertical current sheet above the hot flare loops and below 
the coronal mass ej ection associated with eruptive flares. These and related observations are 
discussed further in lFletcheretaT]d201lh . 
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These recent observations strongly support the "standard" model of eruptive solar flares, 
in which the hot flare loops build up below a vertical current sheet where inflowing magnetic 
fields reconnect an d a magnetic flux rope forms above the cur rent sheet to become a coronal 
mass ejection (see lFletcher et alj|201 1| ; Izharkova et alj|201lh . The rate of electron acceler- 
ation has been observed to be correlated with the rate at which magnetic flux is swept up 
by th e expanding footpoints of flare l oops and with the rate of looptop expansion dOiu et all 
|2004 ISui et al.ll2004l ; lHolmanll200l) . indicating that the electron acceleration rate is cor- 
related with the rate of magnetic reconnection. On the other hand, the observations also 
indicate that the rate of electron acceleration in the impulsive phase of flares is greatest 
before a large-scale current sheet or soft X-ray cusp is observed dSui et al]l200ot) . 

Initially, when the electron acceleration rate is highest, the current sheet may be short 
and assoc iated with slow-m ode shock waves, as in Petschek reconnection. Fast reconnection 
jets (e.g., Wa ng et al.l2007t) can stream upward and downward from the current sheet, possi- 
bly ending in fast-mode shock waves where they collide with slower magnetized plasma at 
the flare loop tops and the lower boundary of the magnetic flux rope (termination shocks). 
The pair of above-the-looptop X-ray sources may be associated with these fast-mode shock 
waves. We have described possible evidence for these shock waves from radio observations 
in SectionlTTl 

The most difficult task is determining the dominant acceleration mechanism or mech- 
anisms responsible for the energetic particles. The region above the flare loops contains or 
can contain quasi-DC electric fields, plasma turbulence, slow- and fast-mode shock waves, 
and collapsing magnetic traps, allowing for almost any acceleration mechanism imaginable. 
Th e problem is as mu ch one of ruling out mechanisms as of finding mechanisms t hat w ork 
Acceleration mechanisms are addressed in lZharkova et al.l d201 lh . 

In Section|9]we addressed the soft-hard-soft evolution of flare X-ray spectra. This spec- 
tral evolution could occur during the propagation of the electrons from the acceleration 
region to the thick-target footpoints. Return current losses, with their dependence on the 
electron beam flux (Section |5j, for example, could be responsible for this evolution. How- 
ever, the observation that above-the-looptop sources also show this spectral evolution (Sec- 
tion ll0.3| > indicates that it is a property of the acceleration process rather than electron beam 
propagation. We saw in Section 19.21 that the soft-hard-soft behavior can be reproduced in 
the acceleration region if the acceleration or trapping efficiency first increases and then de- 
creases. 

Flares displaying soft-hard-harder spectral evolution are of special inter est, because they 
have been shown to be associated with hi gh-energy proton events in space (Kipl ingerj|l995 ; 
ISaldanha et alj|2008l ; lGravson et al.ll2009h . What is the connection between the appearance 
of energetic protons in space and X-ray spectral hardening late in flares? The answer to this 
question is important to both space weather prediction and understanding particle accelera- 
tion in flares. 



12.3 Implications of current results for future flare studies in hard X-rays 

What characteristics should a next-generation hard X-ray telescope have to make substantial 
progress in understanding electron propagation and acceleration in flares? The advances 
made with RHESSI have depended on its high-resolution count spectra that could generally 
be convolved with the detector response to obtain reliable photon flux spectra. These have 
been the first observations to allow detailed information about the evolution of accelerated 
electrons and associated hot flare plasma to be deduced for many flares. Equally important 
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has been the ability to produce hard X-ray images in energy bands determined by the user 
during the data analysis process. This imaging capability has been critical to determining the 
origin of the X-ray emission at a given photon energy and in obtaining spectra for individual 
imaged source regions. These high-resolution imaging spectroscopy capabilities will remain 
important for continued progress. 

RHESSI's X-ray imaging capability has allowed a clear spatial separation to be made for 
many flares between footpoint sources with non-thermal spectra at higher energies and loop- 
top sources with thermal spectra at lower energies. However, in the energy range of overlap 
between ~10 keV and 50 keV, where both types of sources may coexist, it is often difficult 
to distinguish weaker coronal sources (both thermal and non-thermal) in the presence of the 
stronger footpoint sources. This is because of the limited dynamic range of < 100: 1 (and sig- 
nificantly less for weaker events) that is possible in any one image made from RHESSI data. 
This is a consequence of the particular form of the Fourier-transform imaging technique that 
is used. Thus, in most flares the usually intense footpoints mask the much weaker coronal 
hard X-ray source s that can sometimes be seen in over-the-limb flares when the footpoints 
are occulted (e.g., iKrucker & LmteOQSjl In fact, these coro nal hard X-ray sources can ex- 
tend to high energies (up to ~800 keV. lKrucker et alJl2008bh and seem to be non-thermal in 
origin, thus making them of great interest in locating and understanding the particle acceler- 
ation process. It is important to study these non-thermal coronal sources in comparison with 
the footpoint sources, something that is currently not possible with RHESSI's limited dy- 
namic range except in the few cases with exceptionally strong coronal emission (see Section 
10). In addition, again because of the RHESSI dynamic range and sensitivity limitations, it 
has not generally been possible to observe the thin-target bremsstrahlung emission from the 
corona that must be present from the electrons streaming down the legs of magnetic loops 
and also from electrons streaming out from the Sun and producing type-Ill bursts (see, how- 
ever, [^rackeTeyiniooi^ISai^ For all of these reasons, a significantly 
greater dynamic range will be an important goal for future advanced solar hard X-ray in- 
struments. 

Flares at the solar limb for which the hard X-ray footpoints are occulted by the disk 
provide an important way of observing coronal hard X-ray sources, but these flares do not 
allow a comparison to be made between the coronal emission and the thick-target foot- 
point emission. A possible substitute for a high-dynamic-range instrument is hard X-ray 
observations from two or more spacecraft. Under the right conditions, one spacecraft can 
observe all the flare emission while the other observes only the coronal emission, with the 
footpoint emission occulted by the solar disk. Multi-spacecraft observations would also be 
important for deducing the directivity o f the flare emission (especially in conjunction with 
X-ray polarization measurements - see iKontar et al.ll201lh and 3-D source structure. This 
multi-spacecraft approach, however, limits the number of flares for which the coronal and 
footpoint emissions can be compared. 

Hard X-ray timing studies have provided valuable information about electron propa- 
gation and the location of the acceleration region (Section [SJ. Since the time of flight of 
energetic electrons from a coronal acceleration region to the thick-target loop footpoint is 
typically ~ 10-100 milliseconds, the photon count rate must be high enough to distinguish 
differences in flux on these time scales. Time-of-flight studies have not been successful with 
RHESSI, because of its relatively low collecting area and, therefore, count rate. An instru- 
ment with the collecting area and pulse-pileup avoidance of CG7?0/BATSE, and the imaging 
and spectral resolution of RHESSI, would provide a new generation of studies on the char- 
acteristic time scales of propagation for the hard-X-ray-emitting electrons accelerated in 
flares. Alternatively, smaller instruments sent closer to the Sun on, for example, Solar Or- 
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biter or Solar Probe could achieve the required sensitivity. Flare studies on these time scales 
would provide important insights into the physical processes that impact the acceleration 
and propagation of energetic electrons in flares. 
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